Criticality in quantum triangular antiferromagnets via fermionized vortices 
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We revisit two-dimensional frustrated quantum magnetism from a new perspective, with the aim of exploring 
new critical points and critical phases. We study easy-plane triangular antiferromagnets using a dual vortex 
approach, fermionizing the vortices with a Chern-Simons field. Herein we develop this technique for integer- 
spin systems which generically exhibit a simple paramagnetic phase as well as magnetically-ordered phases with 
coplanar and collinear spin order. Within the fermionized-vortex approach, we derive a low-energy effective 
theory containing Dirac fermions with two flavors minimally coupled to a (7(1) and a Chern-Simons gauge 
field. At criticality we argue that the Chern-Simons gauge field can be subsumed into the (7(1) gauge field, and 
up to irrelevant interactions one arrives at quantum electrodynamics in (2+1) dimensions (QED3). Moreover, 
we conjecture that critical QED3 with full SU (2) flavor symmetry describes the 0(4) multicritical point of the 
spin model where the paramagnet and two magnetically-ordered phases merge. The remarkable implication is 
that QED3 with flavor SU(2) symmetry is dual to ordinary critical $ 4 field theory with 0(4) symmetry. This 
leads to a number of unexpected, verifiable predictions for QED3. A connection of our fermionized-vortex 
approach with the dipole interpretation of the v = 1/2 fractional quantum Hall state is also demonstrated. The 
approach introduced in this paper will be applied to spin- 1/2 systems in a forthcoming publication. 

PACS numbers: 



I. INTRODUCTION 

Geometrically frustrated two-dimensional (2D) quantum 
spin systems provide promising settings for the realization of 
exotic quantum phenomena due to the frustration-enhanced 
role of quantum fluctuations. The quasi-2D spin- 1/2 trian- 
gular antiferromagnet CS2CUCI4 and the organic Mott in- 
sulator k-(ET) 2 Cu2(CN)3 provide two noteworthy examples 
that have generated much interest. While CS2CUCI4 exhibits 
"cycloid" magnetic order at low temperatures, evidence of 
a spin-liquid phase with deconfined spinon excitations was 
found at higher temperatures. 1,2 Moreover, application of an 
in-plane magnetic field suppressed magnetic order down to 
the lowest temperatures achieved, which might indicate a sta- 
bilization of the spin-liquid phase. 1 In the organic material 
k-(ET) 2 Cu2(CN)3, NMR measurements show no signs of 
magnetic ordering, suggesting a spin-liquid phase as a can- 
didate ground stated 

Understanding critical spin-liquid states in 2D frustrated 
spin- 1/2 systems is a central issue that has been the focus of 
significant theoretical attention in response to the above exper- 
iments. Several authors have attempted to access spin-liquid 
phases in these systems using slave-fermion and slave-boson 
mean field approaches, 4 ' 5 ' 6 and recent work 7 has studied uni- 
versal properties near the quantum phase transition between 
a spiral magnetic state and one such spin liquid with gapped 
bosonic spinons. Spin-liquid states were also recently exam- 
ined using slave-particle mean field approaches in the trian- 
gular lattice Hubbard model near the Mott transition, which is 
argued to be relevant for k-(ET) 2 Cu2(CN)3. 8 ' 9 

The purpose of this paper is to introduce a new dual vortex 
approach to frustrated quantum spin systems that allows one 
to gain insight into critical spin-liquid phases coming from 
the easy-plane regime. Our goal will be to illustrate how one 
can formulate and analyze an effective low-energy dual the- 



ory describing criticality in frustrated quantum systems using 
fermionized vortices. As we will see, an advantage of this ap- 
proach is that critical points and critical phases arise rather 
naturally from this perspective. 

In the present paper we develop the fermionized-vortex ap- 
proach by studying easy-plane integer-spin triangular anti- 
ferromagnets, postponing an analysis of the more challeng- 
ing spin- 1/2 systems for a future publication. In the integer- 
spin case, a direct Euclidean formulation of the quantum spin 
model is free of Berry phases, and therefore equivalent to a 
classical system of stacked triangular antiferromagnets. Thus, 
a fairly complete analysis is accessible in the direct spin lan- 
guage. Generically, such systems exhibit a paramagnetic 
phase, a magnetically-ordered phase with collinear spin or- 
der, and the well-known coplanar v3 x y/3 state, and much is 
known about the intervening phase transitions. 

Using the direct analysis as a guide, we are led to a num- 
ber of surprising results in the fermionized-vortex formula- 
tion, which we now summarize. We conjecture that criticality 
in the spin model is described by critical quantum electrody- 
namics in (2+1) dimensions (QED3) with two flavors of two- 
component Dirac fermions minimally coupled to a noncom- 
pact gauge field. In particular, we conjecture that this criti- 
cal QED3 theory describes the 0(4) multicritical point of the 
spin model where the paramagnet merges with the coplanar 
and collinear ordered phases. In other words, we propose that 
critical QED3 with SU (2) flavor symmetry is dual to 0(4) 
critical $ 4 theory. Some remarkable predictions for the for- 
mer theory follow from this conjecture, and can in principle 
be used to test its validity. Namely, we predict that single- 
monopole insertions in the QED3 theory, which add discrete 
units of 2ir gauge flux, have the same scaling dimension as the 
ordering field in the 0(4) theory. In fact, the correspondence 
between the two theories can be stated by observing that there 
are two distinct monopole insertions in the QED3, associating 
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a complex field with each monopole, and grouping the two 
complex fields into an 0(4) vector. Furthermore, fermionic 
bilinears ipTip, which are not part of any conserved current, 
and double-strength monopoles in QED3 are predicted to have 
identical scaling dimensions, equal to that of the quadratic 
anisotropy fields at the 0(4) fixed point. It would be quite 
remarkable if the scaling dimensions for these seemingly un- 
related operators indeed merge at the QED3 fixed point. 

Given the length of this paper, it will be worthwhile to now 
provide an overview of the fermionized-vortex approach and 
describe the main steps that lead to our results for the integer- 
spin triangular antiferromagnet. We will first illustrate how 
one obtains a low-energy description of fermionized vortices, 
and then present a more detailed summary of the integer-spin 
case. We will also mention below an application of this ap- 
proach in a different context, namely treating bosons at filling 
factor v = 1, which reproduces the so-called "dipole picture" 
of the compressible quantum Hall fluid. 

We start by highlighting the difficulties associated with 
frustration encountered in a conventional duality approach. 
Frustrated easy-plane spin models are equivalent to systems 
of bosons hopping in a magnetic field inducing frustration, 
and can be readily cast into the dual vortex representation^ A 
crucial feature of the resulting dual theory is that the vortices 
are at finite density due to the underlying frustration. The fi- 
nite vortex density together with the strong vortex interactions 
present a serious obstacle in treating the bosonic vortices, as 
it is not clear how one can obtain a low-energy description 
of such degrees of freedom that describes criticality in the 
spin system. However, we will demonstrate that fermioniz- 
ing the vortices via Chern-Simons flux attachment allows one 
to successfully construct a low-energy dual vortex description 
in both the integer-spin and spin- 1/2 cases. From the view of 
spin systems the principal gain here is the ability to describe 
the latter, as a direct analysis in the spin- 1/2 case is plagued 
by Berry phases. Nevertheless, even in the case of integer-spin 
systems, the dual approach leads to interesting and unexpected 
predictions as mentioned above. 

The formulation of a low-energy dual theory proceeds 
schematically as follows. Let us focus on the easy-plane tri- 
angular lattice antiferromagnet (either integer-spin or spin- 
1/2), which dualizes to a system of half-filled bosonic vortices 
on the honeycomb lattice with "electromagnetic" interactions. 
First, the bosonic vortices are converted into fermions carry- 
ing 2tt fictitious flux via a Jordan- Wigner transformation 11 . 
At this stage, the dual theory describes fermionic vortices 
hopping on the dual lattice, coupled to a Chern-Simons field 
and U(l) gauge field that mediates logarithmic vortex inter- 
actions. Next, one considers a mean-field state in which fluc- 
tuations in the gauge fields are ignored and the flux is taken as 
an average background. Such "flux-smeared" treatments have 
been successful in the composite-particle approach to the frac- 
tional quantum Hall effect. 12 Pursuing a similar approach in 
the context of spin systems, Lopez et alH developed a theory 
of the square-lattice Heisenberg antiferromagnet in terms of 
fermions coupled to a Chern-Simons field. Compared with the 
latter application, we have an advantage working with vortices 
as there is an additional U(l) gauge field which gives rise to 



long-range interactions that strongly suppress vortex density 
fluctuations. This fact will play a central role in our analysis. 
Exchange statistics are consequently expected to play a lesser 
role compared with the case of weakly-interacting fermions. 
Replacing the Chern-Simons flux by its average should then 
be a good starting point for deriving a low-energy theory, and 
the resulting continuum theory should be more tractable. 

It is straightforward to diagonalize the "flux-smeared" 
mean-field hopping Hamiltonian, which at vortex half-filling 
reveals two Dirac nodes in the integer-spin case and four Dirac 
nodes in the spin- 1/2 case. Focusing on excitations in the 
vicinity of these nodes and restoring the gauge fluctuations, 
one finally arrives at a low-energy Dirac theory with either 
two or four species of fermions coupled to the Chern-Simons 
and U(l) gauge fields. In the absence of the Chern-Simons 
field, the theory becomes identical to QED3. 

In the spin- 1/2 case, this dual theory has a very rich struc- 
ture, and describes either a critical point in the spin model or a 
new critical spin-liquid phase. As mentioned above, we defer 
an analysis of spin- 1/2 systems to a subsequent paper, focus- 
ing instead on the integer-spin case. Integer-spin systems pro- 
vide an ideal setting for us to establish the dual approach in- 
troduced here, as much is known about their critical behavior 
from a direct analysis of the spin model. Viewed another way, 
using the known results from a direct approach leads us to the 
remarkable predictions described above for the dual critical 
theory, which are highly nontrivial from the latter perspective. 
We proceed now by providing a more detailed summary of the 
dual analysis for the integer-spin case, beginning with the spin 
model to make the preceding discussion more concrete. 

Throughout the paper, the easy-plane spins are cast in terms 
of quantum rotors. The integer-spin Hamiltonian with nearest- 
neighbor interactions then reads 

H X Y =J^2 C ° S ^ r ~~ ^ r ') + U Y1 U r> (1) 
<rr') r 

where n r is the density, ip r is the conjugate phase, and J,U > 
0. For small U/J [i.e., low temperatures in the corresponding 
(2+l)D classical system], this model realizes the well-known 
V3 x v3 coplanar ordered phase, although deformations of 
the model can give rise to collinear magnetic order as well. 
For large U/J the paramagnetic phase arises. The dual vortex 
Hamiltonian obtained from Eq. Q can be conveniently writ- 
ten 

Hdual = -24 C0S((9 X - X / - Chad) 

(XX') 

+ ]T(VVx - l/2)^ xx ,(iV x , - 1/2) 

xx' 

+ U'^Axafi + J' 4x', (2) 

r (xx') 

with J', U' > 0. Here a XX ' 6 R is a gauge field living on the 
links of the dual honeycomb lattice and e xx ' is the conjugate 
electric field, e 10 ^ is a vortex creation operator, and N x is the 
vortex number operator. (A x a) r in the last line denotes a lat- 
tice curl around each hexagon of the dual lattice. The t v term 
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allows nearest-neighbor vortex hopping, and V xx > encodes the 
logarithmic vortex interactions. This form of the dual Hamil- 
tonian makes manifest the half-filling of the vortices. 

To proceed, the vortices are converted into fermions cou- 
pled to a Chern-Simons gauge field. Because the vortices are 
half-filled, the Chern-Simons flux through each dual plaque- 
tte averages to 2ir, which on the lattice is equivalent to zero 
flux. Hence, in the flux-smeared mean-field state we simply 
set the gauge fields to zero. We emphasize again that this is 
expected to be a good approximation due to the vortex interac- 
tions, which strongly suppress density fluctuations. One then 
has fermions at half-filling hopping on the honeycomb lattice. 
Diagonalizing the resulting hopping Hamiltonian, one finds 
two Dirac nodes. These nodes occur at the same wavevectors 
- namely, the \/3 x \/3 ordering wavevectors - that character- 
ize the low-energy physics in the Landau approach to the spin 
model, which are not apparent in the dual theory of bosonic 
vortices but arise naturally in the fermionized representation. 
A low-energy description can be obtained by focusing only on 
excitations near these nodes, and upon restoring gauge fluctu- 
ations one finds the following Euclidean Lagrangian density, 

1 i 
+ -^{e-^uxduax) 2 + —e^xA^Ax + £ 4 f- (3) 

Here V'fi/i denote the two-component spinors corresponding 
to each node, the flavor index a = R/L is implicitly summed, 
and A is the Chern-Simons field. The usual form for the 
Chern-Simons action above enforces the flux attachment that 
restores the bosonic exchange statistics to the vortices. Fi- 
nally, £ 4 f represents four-fermion terms arising from short- 
range interactions allowed in the microscopic model. 

The phases of the original spin model correspond to mas- 
sive phases for the fermions in this formulation. The para- 
magnetic phase is realized in the dual fermionic language as 
an integer quantum Hall state, with gapped excitations corre- 
sponding to the magnon of the paramagnet. The magnetically- 
ordered states correspond to spontaneous mass generation 
driven by fermion interactions; the \/3 x \/3 spin state 
is obtained as a vortex "charge-ordered" state, while the 
collinearly-ordered states are realized as vortex valence bond 
solids. 

Since the fermionized dual formulation faithfully captures 
the phases of the original spin model, it is natural to ask 
about the phase transitions in this framework. Our explo- 
ration of criticality begins by implementing the following 
trick. First, note that the fermions couple to the combination 
<V = a n + . Rewriting the Lagrangian in terms of this sum 
field, we then integrate out the Chern-Simons field to arrive at 
a Lagrangian of the form 

£ = ^ a 7 M (<9p - ia,j)i) a + -^(e^xd^ax) 2 

+ C if + i-^(V x a) • (V x V x a) . (4) 

Remarkably, the lowest-order term generated by this pro- 
cedure is the last Chern-Simons-like term containing three 



derivatives, which at criticality is irrelevant by power count- 
ing compared to the Maxwell term. It is therefore tempting to 
drop this term altogether, which is tantamount to discarding 
the bosonic exchange statistics and replacing the bosonic vor- 
tices by fermions without any statistical flux attachment, or 
more precisely a complete 'screening' of the statistical gauge 
field. Such a replacement leads to an additional symmetry cor- 
responding to the naive time reversal for fermions without any 
Chern-Simons field, a symmetry which the higher-derivative 
term violates. The unimportance of exchange statistics and 
concomitant emergence of this additional time-reversal sym- 
metry at criticality seems plausible on physical grounds since 
the logarithmic vortex interactions strongly suppress density 
fluctuations. Indeed, the trick performed above is only possi- 
ble because of the additional U(l) gauge field that mediates 
these interactions. We postulate that it is indeed legitimate 
at criticality to drop this higher-derivative term, and that the 
resulting QED3 theory, 

C-QED3 = tpaJ^idfj, - ia^ipa + -^(e^xd^ax) 2 + £ 4 /, (5) 

contains a description of criticality in the spin model. Fur- 
thermore, we boldly conjecture that the critical QED3 theory 
with the full SU(2) symmetry describes the 0(4) multicrit- 
ical point of the spin model where the paramagnet and two 
magnetically-ordered phases meet. 

We support the above conjecture by establishing the 
symmetry-equivalence of operators in the spin model with op- 
erators in QED3. Specifically, we identify the two complex 
order parameter fields of the spin model with the two lead- 
ing monopole insertions of QED3, which add discrete units 
of 2tt gauge flux. Moreover, we establish the equivalence of 
the various spin bilinears with fermionic bilinears and double- 
strength monopole operators of QED3. At the multicritical 
point, all of the spin bilinears (except the energy) must have 
the same scaling dimension due to the emergent global 0(4) 
symmetry. This leads us to another bold prediction: specific 
fermionic bilinears and double-strength monopole insertions 
in QED3 with two fermion species have the same scaling di- 
mension. We again remark that this would be quite surprising, 
as from the perspective of the dual theory there is a priori no 
apparent relation between these operators. 

The above predictions are in principle verifiable in lattice 
QED3 simulations. 14 However, this likely requires fine-tuning 
of the lattice model to criticality, since we expect that the 
system generically ends up in a phase with spontaneously- 
generated mass driven by the fermion interactions that become 
relevant for the low number of fermion flavors in our theory. 
In particular, based on the correspondence with the 0(4) fixed 
point, we expect two four-fermion interactions to be relevant 
at the QED3 fixed point (one strongly relevant and one weakly 
relevant). We note that we have performed a leading-order 
"large-iV" analysis to assess the stability of the QED3 fixed 
point in the limit of a large number of fermion flavors, where 
gauge fluctuations are suppressed. This analysis suggests that 
only one interaction in the two-flavor theory is relevant; how- 
ever, higher-order corrections could certainly dominate and 
lead to two relevant interactions as predicted above. 
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We now want to give a physical view of the fermionized 
vortices and also point out a connection with the so-called 
dipole interpretation of a Fermi-liquid-like state of bosons at 
filling factor v = \^±L This system was studied in connec- 
tion with the Halperin, Lee, and Read (HLR) theory— of a 
compressible fractional quantum Hall state for electrons at 
v = 1/2. In the direct Chern-Simons approach, the two sys- 
tems lead to very similar theories, but many researchers have 
argued the need for improvement over the HLR treatment (for 
a recent review see Ref. [l8l) . A detai led picture for bosons 
at v = 1 is developed in Refs. Il5ll6l working in the lowest 
Landau level. The present approach of fermionizing vortices 
rather than the original bosons provides an alternative and per- 
haps more straightforward route for describing this system 
and obtaining the results of Ref. 16. Our approach is closest in 
spirit to Ref. 19, which also uses duality to represent vortices. 
The fermionized vortices appear to be 'natural' fermionic de- 
grees of freedom for describing the compressible quantum 
Hall state. We also note an early work 20 where the qualita- 
tive possibility of fermionic vortices was raised in a context 
closer to frustrated easy-plane magnets. 

The vortex fermionization procedure is essentially the same 
as for the spin model described above. More details for this 
case can be found in Appendix [A] Dualizing the v = 1 
bosons, one obtains a theory of bosonic vortices at filling fac- 
tor i^duai = 1/v = 1, with mean vortex density equal to the 
density of the original bosons p. Again, unlike the original 
bosons, the vortices interact via a 2D electromagnetic inter- 
action, and fermionization of vortices gives a more controlled 
treatment. We obtain fermions in zero average field, coupled 
to a U(l ) gauge field, a, and a Chern-Simons gauge field, A. 
This theory is very similar to Eq. (0, but with nonrelativistic 
fermions at density p. With this theory, one can, for exam- 
ple, reproduce the physical response properties expected of 
the compressible state — by using an RPA approximation for 
integrating out the fermions. 

It is useful to exhibit the equivalent of Eq. @ in the present 
case, 

C v=l = V f (d T - [ ^~ m] ) ip - ia Q {i>H -p) + ..(6) 

Here m VO rt is a phenomenological vortex mass. As before, 
a = a + A, and the Chern-Simons field has been integrated 
out. The final action has nonrelativistic fermions coupled to 
the gauge field a. As detailed in Appendix |A] in the lowest 
Landau level limit there is no Chern-Simons like term for the 
field a in contrast to Eq. @, and the theory has the same struc- 
ture as in Ref. 

In the present formulation, the fermions are neutral and do 
not couple directly to external fields. However, their dynamics 
dictates that they carry electrical dipole moments oriented per- 
pendicular to their momentum k and of strength proportional 
to |fe|, which can be loosely viewed as a dumbbell formed by 
an original boson and a vortex. This can be seen schemati- 
cally as follows. First we note that the flux of a is equal to the 
difference in the boson and vortex densities: 

V Ao = V Aa + V A A = 2n(p hos - p vort ) . (7) 



Here VAq = Cijdihj denotes a 2d scalar curl; we also define 
(/\v)i — CijVj for a 2d vector v. Consider now a wavepacket 
moving with momentum k in the y direction, and imagine a fi- 
nite extent in the x direction. In the region of the wavepacket, 
we have a y — \k\, while a y is zero to the left and to the right 
of the wavepacket. Then V Aa ~ d x a y is positive near the left 
edge and negative near the right edge, i.e. the boson density 
is higher than the vortex density on the left side of the moving 
fermion and vice versa on the right side. A more formal cal- 
culation is presented in Appendix 1X1 and gives the following 
dipolar expression for the original boson density in terms of 
the fermionic fields: 

Sp bos (q) =iq- j ^{- q/2 ^k+ q /2 ■ (8) 

Thus, by fermionizing the vortices we indeed obtain a de- 
scription of the compressible v — 1 state very similar to the 
dipole picture, and this is achieved in a rather simple way 
without using the lowest Landau level projection. The suc- 
cess of our approach in this case makes us more confident in 
the application to frustrated spin systems. 

The remainder of the paper focuses on the integer-spin tri- 
angular antiferromagnet and is arranged as follows. In Sec. 
HHwe review the phases and critical behavior that arise from 
a direct analysis of the spin model. The duality mapping is 
introduced in Sec. |ni] In Sec. lIVI the fermionized-vortex for- 
mulation is developed. A low-energy continuum theory is ob- 
tained here, and we also discuss how to recover the phases 
of the spin model from the fermionic picture. An analysis of 
the QED3 theory conjectured to describe criticality in the spin 
model is carried out in Sec. [V] We conclude with a summary 
and discussion in Sec. lVIl 

II. DIRECT LANDAU THEORY 

A. The Model 

In this paper we focus on quantum triangular XY antiferro- 
magnets with integer spin, described by the Hamiltonian 

F = J/2 ^ (S+S- + S-S+) + Uj2(Sr) 2 , (9) 

<rr') r 

with J, U > 0, S± = S?±iS$, and [S?, Sf,] = i6 rr >e a07 S? . 
Rather than working with spin models, it will be particularly 
convenient to instead model such easy-plane systems with ro- 
tor variables, introducing an integer-valued field n r and a 2tt- 
periodic phase variable ip T that satisfy [ip r , n T >] = i5 rr >. Upon 
making the identification, 5+ — > e liPr and S* — > n r , the ap- 
propriate rotor XY Hamiltonian reads, 

(rr'> r 

Equation dlOt can equivalently be written 

HXY = -jJ2 C ° S ^ r - ^' + A ™' ) + U J2 U r> (U) 
<rr') r 
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where we have introduced new phase variables <fi r and a static 
vector potential A® T , that induces tt flux per triangular pla- 
quette. In this form, the nearest-neighbor couplings can be 
viewed as ferromagnetic, with frustration arising instead from 
the flux through each plaquette. The model can alternatively 
be viewed as bosons at integer-filling hopping on the trian- 
gular lattice in an external field giving rise to one-half flux 
quantum per triangular plaquette. 

A closely related model also described by Eq. (II It is the 
fully-frustrated XY model on a square lattice, which has ap- 
plications to rectangular Josephson junction arrays. The tech- 
niques applied to the triangular antiferromagnets in this paper 
can be extended in a straightforward manner to this case as 
well. 

Within a Euclidean path integral representation, the 2D 
quantum model becomes identical to a vertical stacking of 
classical triangular XY antiferromagnets with ferromagnetic 
XY coupling between successive planes, imaginary time giv- 
ing one extra classical dimension. This classical model has 
been extensively studied, and much is known about the phases 
and intervening phase transitions. We will be particularly in- 
terested in the phase transitions, which are quantum transi- 
tions in the spin model. In this section we will briefly review 
the known results obtained from a Landau-Ginzburg-Wilson 
analysis and from Monte Carlo studies. This summary of the 
direct approach to these models will later serve as a useful 
counterpart to our dual approach. 



B. Low-energy Landau Theory 

A low-energy effective theory for the triangular XY an- 
tiferromagnet can be obtained by introducing a Hubbard- 
Stratonovich transformation that replaces e lVi by a "soft spin" 
with the same correlations. This leads to a Euclidean ac- 
tion of the form 

S = I []r(|^$ r | 2 +m|<i> r | 2 + U |<i> r | 4 ) 



<rr') 



(12) 



with J > 0. [Note that we used the angle variables appear- 
ing in Eq. dlOl to implement this transformation.] Diagonal - 
izing the kinetic part of the action in momentum space, one 
finds two global minima at wavevectors ±Q, where Q = 
(47r/3, 0). The low-energy physics can be captured by ex- 
panding the fields around these wavevectors by writing 



$(r) 



*fl(r) 



<Mr), 



(13) 



where Qr/l are complex fields assumed to be slowly-varying 
on the lattice scale. 

To explore the phases of the system within a Landau- 
Ginzburg-Wilson framework, one needs to construct an effec- 
tive action for these fields that respects the microscopic sym- 
metries of the underlying lattice model. In particular, the ac- 
tion must preserve the following discrete lattice symmetries: 



translation by triangular lattice vectors Sr (Ts r ), ^-reflection 
[1Z X : (x,y) — > (— x, y)], and 7r/3 rotation about a lattice 
site (R n / 3 ). Additionally, the action must be invariant under 
a number of internal symmetries, specifically a global U(l) 
symmetry, (p r — > <p r + a; a "particle-hole" or "charge" conju- 
gation (C) symmetry, 

C : n r —> -n T , tp r -> -tp r ; (14) 

and spin time reversal (T) which sends — ► —S", 

T : n r — ► — n r , tp T — > (fir + 7r, i — > —i. (15) 

Under these symmetry operations, the continuum fields ^r/l 
transform according to 



Tst 

U(l) 

c 

T 



R/L 
§R ®L 



e ±iQ - Sr $> R /L, 



R/L 
R/L 
R/L 



R/Li 



(16) 



L/R> 



In the last line above, we followed time reversal by a U(l) 
transformation to remove an overall minus sign. 

These symmetry constraints lead to a continuum action of 
the form 

S *ff = / { E (lW + r|*a| a ) 

•>t,t I „ it 



a=R/L 

«4(|*H| 2 + |^i| a ) 2 + «4|*il| 2 |^£| S 



^Wr^l? + C.C.)] 



(17) 



The v§ term has been included since it is the lowest-order term 
that reduces the symmetry down to the global U(l) of the 
microscopic model. Note that when V4 = v$ = the action 
acquires an 0(4) symmetry. On the other hand, when v§ = 
and V4 = —2u4 the action reduces to two decoupled U(l) 
models. 

We now proceed to explore the phases that arise from 
Eq. d!7l > within mean-field theory, assuming U4 > but al- 
lowing r, V4, and vq to be either positive or negative. 



C. Phases 

1. The simplest phase arises when r > so that ($ a ) = 0. 
This corresponds to the quantum paramagnet with (S~) = 
at each site. 

When r < 0, magnetic order develops, with the character 
depending on the sign of V4: 

2. The case r < 0, V4 > favors a state with ($r) 7^ and 

= (or vice versa). Taking = 1, the order param- 
eter in terms of the original spin variables is (S~) = e lC *' T . 
This corresponds to the y3 x V3 state with 120° coplanar 
order, whose spin configuration is illustrated in Fig. [2a). 
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3. When r, v A < 0, states with ($r/l) ~ e i0R / L ^ 
emerge. These states are characterized by collinear spin order, 
which can be obtained from 

(S~) ~e ie+/2 cos(Q-r + 6L/2), (18) 

where 9± = Or . ± &l. Two types of collinear spin states 
can occur depending on the phase difference 0_, which is 
determined by the sign of v%. (i) For vq > 0, the action is 
minimized with a phase difference of 9- = 2mr/3, where 
n is an integer. The resulting spin order on the three sublat- 
tices of the triangular lattice can be written schematically as 
(1, — |, — i). This can be viewed as the dice-lattice collinear 
antiferromagnetic state depicted in Fig. \l[b), where "happy" 
nearest-neighbor bonds with spins aligned antiparallel form a 
dice lattice as shown by the solid lines. There are three distinct 
such ground states (in addition to broken global U(l) sym- 
metry) corresponding to the three inequivalent values of . 
(ii) When v$ < 0, a phase difference of = (2n + l)7r/3 
is preferred. The spin order on the three sublattices is then 
(1, — 1, 0), which can be viewed as hexagonal collinear anti- 
ferromagnetic order as illustrated in Fig. Qc). In this state, 
spins on one of the three hexagonal subnetworks of the tri- 
angular lattice order antiferromagnetically, while spins on the 
remaining sites fluctuate around zero average. Here also there 
are three distinct ground states arising from the inequivalent 
values of . 

One can also distinguish between the coplanar and collinear 
magnetic orders by considering chirality and bond energy 
wave operators. Spin chirality k p at each triangular plaque- 
tte p can be defined by 

K P ~ 1C sin (Vr' - ¥>r), (19) 
(rr'J gp 

where the sum is over the links of the plaquette oriented coun- 
terclockwise. A local chirality operator is obtained by sum- 
ming k p over local "up" triangles. A bond energy wave oper- 
ator may be defined as 

B rr , ~ e 4Q ( r+r ') cos((^ r - yv) . (20) 
In terms of the continuum fields $>r/l> we identify 

k ~ \<S>r\ 2 - |$l| 2 =K z , 
B ~ <S>* R $ L = (K x + iK y )/2. (21) 

Here we have introduced a three-component vector field, 

K = $;r Qb $ b , (22) 

with t = (t x ,t v ,t z ) a vector of Pauli matrices. In the copla- 
nar phase, the vector K points in the ±K Z direction, while in 
the collinear phase it lies in the (K x , K y ) plane. The degener- 
acy in the plane is lifted by the vq term, leaving three possible 
ordering directions 6>_ for each collinear phase as discussed 
above. 

Microscopically, the coplanar, collinear, and paramagnetic 
phases can be realized in the classical stacked triangular 



(a) 

Coplanar Spin State 




FIG. 1 : (a) Coplanar spin configuration in the v3 x v3 ordered state, 
(b) Dice-lattice collinear antiferromagnetic state. Bold lines denote 
satisfied bonds, which form a dice lattice, while the remaining bonds 
are frustrated. | {<S'~ ) | is roughly twice as large on sites with long ar- 
rows than on sites with short arrows, (c) Hexagonal collinear antifer- 
romagnetic state. Circles represent fluctuating spins with (S~) — 0. 
Satisfied bonds denoted by bold lines form a hexagonal lattice. 

XY antiferromagnet with an additional "boson pair hopping" 
term, 

Hd = COs[ip ryZ + 1 - (f T>z ] 

r.z 

+ J ^ COS^r,* - Vr'.z] 
(rr') ,z 

- J' ^ COs[2(^r,z - <Pr',z)], (23) 
<rr') : z 

where J' > and the z coordinate labels the different tri- 
angular lattice planes. The y/3 x v3 coplanar phase is the 
ground state when J' <C J. As J' increases, collinear anti- 
ferromagnet order eventually becomes more energetically fa- 
vorable. Of course, at high enough temperature we obtain the 
paramagnetic phase. A sketch of the phase diagram contain- 
ing the three phases is shown in Fig.|2ja). 



D. Critical Behavior 

Significant numerical and analytical effort has been di- 
rected towards understanding the phase ttansitions occurring 
in classical stacked triangular XY antiferromagn ets. Fo r a re- 
cent review, the reader is referred to Refs. 12 Il22l23l24l Here, 
we shall only highlight the known results that are pertinent to 
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(a) (b) 




r 



FIG. 2: (a) Sketch of the phase diagram in the vicinity of the 
coplanar-collinear transition, (b) Putative renormalization group 
flows in the (r, V4) plane assuming a continuous transition from both 
the coplanar and collinear states to the paramagnet. 



the present work. A brief discussion of our own exploratory 
Monte Carlo study of the classical model Eq. d23i that real- 
izes both the coplanar and collinear (hexagonal) spin-ordered 
phases will also be given. 

In mean field theory, the paramagnet to magnetically- 
ordered transitions are continuous, while the transition from 
the coplanar to the collinear state is first order. 

The transition between the paramagnet and the collinear 
state also appears to be continuous in our Monte Carlo study 
of the model Eq. (1231 near the multicritical point where the 
three phases meet. Moreover, such a continuous transition is 
expected to be governed by a /7(1) x E7(l) fixed point with 
v\ = —2^4 and Wg = 0. Small deviations from the condition 
i<4 = — 2^4 give rise to energy-energy coupling of the two 
U(l) models, which is irrelevant at the decoupled fixed point 
in three dimensions. 22,25 It should be noted that a 4 — e expan- 
sion in the V4 < case predicts a continuous transition charac- 
terized by a different stable fixed point, while the U(l) x U(l) 
fixed point is found to be unstable in that approachi 2 ^ 2 ^ Thus 
in this case the 4 — e expansion does not capture the critical 
physics in three dimensions. 

The nature of the transition between the paramagnetic and 
coplanar phases is controversial! 2 ^ 3 * 2 ^ Starting with Kawa- 
mura, several Monte Carlo simulations on the stacked tri- 
angular XY antiferromagnet conclude that the transition is 
continuous, 2 ^ 2 ^ 2 ^^ which also appears to be the case in our 
Monte Carlo simulations in the vicinity of the multicritical 
point. However, recent simulations on larger lattices claim a 
very weak first-order transition, 31 ' 32 and simulations of modi- 
fied XY systems expected to be in the same universality class 
have also seen a first-order transitions 2 ^ On the other hand, 
Ref. |3 claims a continuous transition in the vicinity of the 
0(4) fixed point in a model which is a lattice discretization of 
the action Eq. ( fTTl . 

Analytical results for the paramagnet-to-coplanar transi- 
tion are conflicting as well. For instance, recent six-loop 
perturbative renormalization group calculations in d = 3 
dimensions 34,35 find a stable "chiral" fixed point as conjec- 
tured by Kawamura 26 . On the other hand, "non-perturbative" 
renormalization group techniques predict a weak first-order 
transition. 23 The 4 — e expansion in the XY case (n = 2 case 



TABLE I: Transformation properties of the bosonic fields and bilin- 
ears under the microscopic symmetries. Here K± = K x ± iK y and 
I± = I x ± il y . For brevity, the lattice coordinates which should 
transform appropriately under the lattice symmetries have been sup- 
pressed. We also show field transformations under the modified re- 
flection TZx = IZxCT, which will be useful later. 







Rtt/3>K x 


C 


T 


Kx 


<$>R/L -<• 




&L/R 


®l/R 


®*R/L 




K x 


K z 


-K t 


-K z 


K z 


K z 


K ± -> 


e =F ^ <4dr if± 




K± 


K T 


K ± 


/. -» 


h 


h 


-11 


7; 


-h 


I± 


e ±M.U-*r J± 


-h 


7* 


4 


-I± 



in Ref. 26) finds runaway flows for V4 > 0, which would be 
interpreted as predicting a first-order transition. 

Finally, the multicritical point where the three phases meet 
is described by an 0(4) fixed point with v\ = Vq = 0; this 
fixed point is unstable towards introducing a small nonzero V4 
term. 

Figure|3b) contains a sketch of the renormalization group 
flows in the (r, V4) plane expected if the transition from the 
paramagnet to the coplanar state is indeed continuous and 
governed by the corresponding fixed point, labeled "C". (In 
the figure, the U4 axis is perpendicular to the plane, and we 
are considering fixed points that are stable in this direction.) 
Fixed point C would disappear if this transition is first-order, 
leading to unstable flows toward V4 — * 00. We note the pos- 
sibility of observing crossover behavior controlled by the un- 
stable 0(4) fixed point shown in Fig.|2jb) by fine-tuning both 
r and V4. This can be explored with Monte Carlo simulations 
of H c i by observing scaling in the vicinity of the multicritical 
point where the paramagnet, coplanar, and collinear phases 
meet in Fig.[2a). 

To examine the 0(4) symmetry at the multicritical point, it 
is convenient to express $b,/l as 

$r =Xi + iX2, ®l = Xa + «X4, (24) 

where \j arc real fields. One expects that \ = 
(Xi> Xii X3j Xi) should transform as an 0(4) vector, and the 
ten independent bilinears XiXj can be decomposed into an 
0(4) scalar, 

K a = $:$ a = x 2 , (25) 

and a traceless, symmetric 4x4 matrix that transforms as a 
second rank tensor representation of 0(4). It will be useful 
later to observe that these nine bilinears can be arranged into 
the SU (2) vector K defined in Eq. \22l . together with two 
other SU(2) vectors: 

/ = $ a [xr y ] ah $ b , (26) 

and its Hermitian conjugate, /*. Under the global U(l) sym- 
metry one has K — > K and I — > e 2la I. For future com- 
parisons the transformation properties of these nine bilinears 
under the remaining microscopic symmetries are shown in Ta- 
ble|I] The 0(4) scalar Kq defined in Eq. j25t is invariant un- 
der all of the above microscopic symmetries. 



8 




FIG. 3: Triangular lattice (dashed lines) and its dual honeycomb lat- 
tice (solid lines). The shaded region indicates the two-site unit cell 
taken for the honeycomb, with the bottom site labeled "1" and the 
top site "2". 

One can furnish explicit microscopic realizations of the bi- 
linears I and I* , just as we did in Eq. (12 1 1 for the bilinears 
K. In our Monte Carlo study of the model Eq. i23\ . we moni- 
tored the apparent scaling dimensions of the nine bilinears and 
find that these indeed merge upon approaching the multicrit- 
ical point, consistent with the expectation for the 0(4) fixed 
point. 

m. DUAL VORTEX THEORY 
A. Duality Mapping 

In this section we introduce the dual formalism that will be 
used throughout the remainder of the paper. At this point, it 
will prove convenient to work with the Hamiltonian as written 
in Eq. fH^, where frustration arises from the flux induced by 
the vector potential A°. r , . We will implement the standard XY 
duality, 10 expressing the Hamiltonian in terms of gauge fields 
living on the links of the dual lattice, which in the case of the 
triangular lattice is the honeycomb (see Fig. [3). 

As a first step, we define oriented gauge fields e xx > E 
[—1/2, 1/2) and a xx < E 2irZ living on the links of the hon- 
eycomb lattice, where x and x' label nearest-neighbor dual 
lattice sites. (We reserve "x" and "r" for sites of the honey- 
comb and triangular lattices, respectively.) The dual "electric 
field" e xx / and "vector potential" a xx > satisfy the commuta- 
tion relation [e XX ', a xx <] = i on the same link and commute 
on different links. 

The boson number and phase are related to the dual fields 
as follows: 

n r = -!-(A x a) r = -!- Wa xx , (27) 

Z7T Z7T z — ' 

(xx')Sp r 

oo 

g#r = JJ. e 27T«W . (28) 

r 

On the right side of Eq. ( I27l i. the sum is over the 
counterclockwise-oriented links of the dual plaquette p r en- 



closing site r of the triangular lattice. The product in Eq. (1281 
is carried out over an arbitrary string running along the trian- 
gular lattice from site r to spatial infinity. A factor of e 2mexx ' 
appears in the product for every link on the dual lattice bi- 
sected by the string, with x on the "right" side of the string 
and x' on the "left". 

To ensure the uniqueness of the phases (j> r obtained from 
Eq. ( l28l . the dual Hilbert space is constrained so that the op- 
erator 

7V X = (A.g) x = ^g xx , (29) 

has integer eigenvalues for all x. On the right-side of Eq. 
J29b . the sum is over the three nearest-neighbor sites of x. 
The meaning of iV x can be understood by using Eq. J28t > to 
show that 

N * = ^ £ for* - (30) 
(rr')ep x 

where the sum is over the counterclockwise-oriented links of 
the triangular plaquette p x enclosing site x and the phase dif- 
ference between adjacent sites is taken to lie in the interval 
[— 7r, 7r). The eigenvalues of this operator therefore encode 
vorticity, i.e., the winding number of the phase around the tri- 
angular plaquette enclosing x. 

To complete the duality transformation, define a static elec- 
tric field e xx , related to the vector potential A°. r , of Eq. (flTt 
by 

e o , = 1 ao ft,) 

l 'xx' 27T ' 

where e xx , and A^ r , live on intersecting links of the honey- 
comb and triangular lattices. Since A® r i generates it flux per 
plaquette of the triangular lattice, we require (A • e°) x to be 
half-integer valued. A convenient gauge choice we will now 
employ is _4" r , = ±7r/3 on each nearest-neighbor link, di- 
rected clockwise around "up" triangles and counterclockwise 
around "down" triangles. In this gauge, e xx , = 1/6 along the 
arrows shown in Fig.[3]so that (A • e°) x = +1/2 on sublattice 
1 while (A • e°) x = -1/2 on sublattice 2. 

Upon defining e xx / = e xx / — e xx / we thereby arrive at a 
simple form for the dual Hamiltonian, 

H = - J ^2 cos(27re xx /) + [/'^](Ax af r , (32) 

<xx') r 

with U' = U/(2ir) 2 . This Hamiltonian is supplemented by 
the constraint that, 

(A • e) x = N X - 1/2, (33) 

where we have found it convenient to shift the integer field 
N x — ► N x — 1 for all x on sublattice 2. 

As it stands, this Hamiltonian is difficult to work with due 
to the integer constraint on the field a xx > / (2ir). It is therefore 
highly desirable to "soften" this constraint, allowing a xx / to 
roam over all real numbers. To be consistent e xx / must also be 
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taken on the real numbers, and it is then legitimate to expand 
the cosine term in Eq. J32l > to obtain, 

H = ^[J'4 X , -2^008^ -0 X / -OxxO] 

(«') 

+ [/'^(Axa)2, (34) 

r 

with J' ps 2ir 2 J. Here we have added a cosine term acting on 
the field a xx < to implement the integer constraint softly, and 
explicitly expressed the longitudinal piece of lattice 
derivative of a "phase" field 6 X S [— 7T, 7r) living on the dual 
lattice sites. The dual Hilbert space satisfies the constraint Eq. 
(I33i with the vortex number operator N x e Z conjugate to 
# x . The operator e * creates a vortex at site x and satisfies 
the commutation relations [N x , e l6x '} = 5 xx >e l6x . The t v term 
above thus allows nearest-neighbor vortex hopping in addition 
to energetically favoring integer values for the density n T . 

Generically, we should also allow short-range vortex inter- 
action of the form 

5](iV x -l/2)y x s x ,(iV x ,-l/2), (35) 

xx' 

in addition to the retarded long-range interaction mediated by 
the gauge fields. For a static configuration of vortices or in 
the instantaneous limit when the "speed of light" ([/'/J') 1 / 2 
is infinite, the latter interaction becomes a 2D lattice Coulomb 
potential varying as V^" 1 ~ — J' In |x — x'|. 

B. Phases in the Dual Vortex Variables 

Physically, the dual Hamiltonian describes vortices hopping 
on the sites of the dual honeycomb lattice, interacting via a 2D 
"electromagnetic" interaction. Most importantly, the vortices 
are at half-filling, which can be directly traced to the frustra- 
tion on the original triangular lattice plaquettes. This high- 
lights the distinction between the triangular lattice XY anti- 
ferromagnet and its unfrustrated counterparts, which have a 
dual vortex theory with an integer mean vortex number. 

To illuminate the challenge in treating the strongly inter- 
acting vortex system when at half-filling, it is instructive to 
consider how the various phases of the original model are 
described in terms of the dual variables. The magnetically- 
ordered states correspond to "insulating" phases of the dual 
vortices. For example, the ordered phase with 120° coplanar 
order of the XY spins corresponds to a "vortex density wave" 
state in which the vortices sit preferentially on one of the two 
sublattices of the dual honeycomb lattice. On the other hand, 
the ordered states with collinear order correspond to "vortex 
valence bond" phases; one signature of these phases is the 
bond energy density order, which can be measured both in 
terms of the vortices and in terms of the original spins. Specif- 
ically, the bond energy wave operator for the vortices takes a 
similar form as in Eq. (I20> . except with the vortex bond energy 
operator inserted, 

B rr , ~ e «*C«+'') cos(0 x - 8 X , - a xx ,) . (36) 



Here, the nearest-neighbor pairs r, r' and x, x' are bridged by 
intersecting links of the triangular and honeycomb lattices. 

In all of these magnetically-ordered phases, since the dual 
vortices are "insulating" and immobile, the "photon" in the 
dual fields a XX ' and e xx / can freely propagate. This "photon" 
mode corresponds to the Goldstone spin-wave mode of the 
original XY spins. On the other hand, the paramagnetic phase 
for the XY spins corresponds to a phase within which the dual 
vortices have condensed, (e l0x ) ^ 0. In this dual vortex "su- 
perfluid" the dual gauge flux is expelled, and the dual gauge 
field a XX ' picks up a Higgs mass; thus, the spectrum is gapped 
in this phase as expected in the spin paramagnet. 

The quantum phase transitions separating the paramagnet 
from the coplanar and collinear ordered spin states correspond 
to "superfluid-insulator" transitions for the dual vortices. But 
all of these vortex "insulators" involve a spontaneous breaking 
of lattice translational symmetries, due to the fact that the vor- 
tices are at half-filling. Due to these complications it is not at 
all clear how one can possibly access the critical properties of 
these transitions from this dual vortex formulation. In the next 
section we demonstrate how this can be achieved, by fermion- 
izing the vortices. Before undertaking this rather tricky busi- 
ness, we first summarize how the vortex fields transform under 
all of the lattice and internal symmetries. 



C. Symmetry Transformations for Vortices 

The transformation properties of the dual vortex and gauge 
fields can be deduced upon inspection of their defining equa- 
tions, Eqs. 1271 and d28i . together with the transformation 
properties of the rotor fields <fi T and n r . For example, lattice 
translations T$ r and tt/3 rotations R n /^ transform the dual 
coordinates (e.g., x — > x + Sr under T$ r ), but do not trans- 
form the dual gauge fields a, e or vortex fields 8, N as shown 
in Table HJ Here and in what follows, our coordinate origin 
is always at a triangular lattice site. Under charge conjugation 
C, however, the fields a,e,8 change sign while N — > 1 — N. 
Similarly, time reversal T changes the sign of the dual fields 
a, 8 but leaves e, N unchanged. Since T is accompanied by 
complex conjugation the vortex creation operator e ld remains 
invariant under time reversal. 

For lattice reflections it is convenient for later developments 
to henceforth consider a modified (antiunitary) transforma- 
tion, 1Z X , which is defined as 1Z X combined with C and T: 

U x = K X CT. (37) 

The dual fields transform under 1Z X in the same way they 
transform under time reversal, with the dual coordinates ad- 
ditionally x-reflected. The complete set of symmetry trans- 
formations for the dual fields is summarized in Table UTI 

Finally, we remark that the global U(l) or XY spin symme- 
try corresponding to e 1 ^ — > e l ^ +a ' is not directly manifest 
in the dual vortex formulation, as can be seen by examining 
Eq. i28l . This global U(l) symmetry reflects the underlying 
conservation of the boson number n r (or the S z -component 
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TABLE II: Transformation properties of the dual fields under the 
microscopic symmetries. For brevity, the lattice coordinates, which 
also transform appropriately under the lattice symmetries, have been 
suppressed on all fields. Here we use a modified reflection TZ X , which 
is a combination of JZ X , C, and T. Both T and 1Z X are accompanied 
with a complex conjugation, so that e.g. e %e is invariant. 



Tfir, Rtt/3 


C 


n x ,T 


a — *• a 


a —* —a 


a —t —a 


e — * e 


e — ► — e 


e — * e 


6^8 


-> -9 


0-» -6» 


TV -> TV 


TV ^ 1 — TV 


TV -> TV 



of spin) and is replaced in the dual formulation by the con- 
servation of the dual flux (A x a) r . Note also that the dual 
theory has a U(l) gauge redundancy (not a physical symme- 
try) given by a xx < -> a xx , + A x - A x , and 9 X -> X + A x , 
where A x G K. 



IV. FERMIONIZATION OF VORTICES 
A. Chern-Simons Flux Attachment 

As described in Sec. IIII Bl an analysis of the phase tran- 
sitions in the dual bosonic vortex theory is highly nontrivial 
because the vortices are at half-filling, and we do not know 
how to formulate a continuum theory of such bosonic degrees 
of freedom. We will sidestep difficulties associated with the 
finite vortex density by fermionizing the vortices, which as we 
will demonstrate enables one to make significant progress in 
this direction. Specifically, we first treat the vortices as hard- 
core bosons, with the identification 

TV X = &t 6 x = o, 1 . (38) 

This is a reasonable approximation due to the repulsive vor- 
tex interactions, and the hard-core condition does not affect 
generic behavior. We then perform a two-dimensional Jordan- 
Wigner transformation!!, 

fct = 4exp[i arg(x,x')iV x ,], (39) 

TVx = 6t&x = <4dx. (40) 

Here, arg(x, x') is an angle formed by the vector x — x' with 
some fixed axis. It is simple to check that d x are fermionic op- 
erators satisfying the anticommutation relations {d Xl d' x ,} = 
Sxx' ■ The vortex hopping part of the Hamiltonian becomes 

-t v J2 [€^- iAx ^ 2 e- ia ^ +h.c.]. (41) 

<xix 2 > 

The Chern-Simons field 

^x lX2 = [arg(x 2 ,x')-arg(x 1 ,x')]A r X ' (42) 

x'^xi,x 2 



TABLE III: Symmetry transformation properties of the fields in the 
fermionized representation. The transformation properties of a, e are 
the same as in Table IIII In the C column, j = 1,2 refers to the 
sublattice index on the honeycomb lattice (site index in the unit cell 
pictured in Fig.[3j- The additional column 7f crm corresponds to the 
naive time reversal for the lattice fermions and is not a symmetry of 
the vortex Hamiltonian. 
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—a 
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—a 


e — ► 
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— e 
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d x -> 
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A — > 


A 


-A 


-A 
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lives on the links of the honeycomb lattice and is completely 
determined by the positions of the particles. In this transfor- 
mation, we have essentially expressed hard-core bosons on the 
lattice as fermions each carrying a fictitious 2tt flux. We used 
the Hamiltonian language in order to facilitate our discussion 
of the discrete symmetries below. 

The Chern-Simons field satisfies 

(A x A) r = y ( N * - V 2 ), (43) 

X£p r 

i.e., the fictitious 2ir flux attached to a fermion can be viewed 
as equally-shared among its three adjacent plaquettes. In the 
above equation, we subtracted 1/2 from each TV X to remove 
an unimportant 2tt flux from each hexagon. This constitutes 
a convenient choice such that the Chern-Simons flux piercing 
the dual lattice vanishes on average since the vortices are half- 
filled. There are further restrictions on the Chern-Simons field 
which we write schematically as A • A = appropriate in the 
continuum limit. 

The complete Hamiltonian describes half-filled fermions 
with nearest-neighbor hopping on the honeycomb lattice, cou- 
pled to the U(l) gauge field a xx > and the Chern-Simons field 
^4xx' ■ We can crudely say that the gauge field a xx / gives rise 
to repulsive logarithmic interactions between the fermions. 

Under the microscopic symmetries, the fields in this 
fermionized representation transform as shown in Table [H3] 
The symmetry properties of the fermion operators under trans- 
lation, rotation, and modified reflection are readily obtained 
by inspecting Eq. ( I39l l (throughout, we ignore a possible 
global U(l) phase). In each case, the transformation of the 
Chern-Simons field is obtained from Eq. J42l >. 

Charge conjugation and time reversal require some expla- 
nation. From Eq. d39> . we have for particle-hole 

C : d x — * d^. exp[i ^ arg(x, x')] = cft x exp[z7 x ] . (44) 

x'^x 

The phase 7 X is constant in time but depends on the lattice 
site. For two nearest-neighbor sites on the honeycomb lattice 
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we have 

7 Xl -7x 2 = arg(xi,x 2 ) - arg(x 2 ,xi) 

+ [ ar S( x i' x ') - ai "g( x 2,x')] 

x'^Xi,X2 

= 7T - 2A xlX2 [iV x , = 1/2] . 

With our convention in Eq. J43> where the Chern-Simons field 
fluctuates around zero, we conclude that e l7x changes sign 
going between nearest neighbors. The resulting e' nx is written 
as (— iy in Table lllll where j = 1,2 refers to the sublattice 
index of site x. 

Finally, time reversal acts as 

T : d x -> d x cxp[-2i ^ arg(x, x')iV x ,]. (45) 

x'^x 

This is a formally exact implementation of the symmetry, but 
the phases multiplying the fermion operators now depend on 
the positions of all the vortices. This nonlocal transformation 
represents a serious difficulty once we derive a low-energy 
continuum theory in the next subsection. In particular, it will 
not be possible to correctly represent the time-reversal trans- 
formation T using only the continuum fields. Since our treat- 
ment below will involve crucial assumptions regarding this 
issue, it is useful to also give a formulation of time-reversal 
symmetry in first-quantized language. Focusing on the par- 
ticle degrees of freedom, the Jordan-Wigner transformation 
expresses the wavefunction for the bosonic vortices as 

*boson(xi,X 2 , . . . ) = JJ e -*«8(x.«) vP fcrm ( Xl ,X 2 , . . . ) . 

i<j 

(46) 

When the Chern-Simons phase factor is not affected by a 
symmetry transformation, the properties of the bosonic and 
fermionic wavefunctions coincide under the transformation. 
However, time reversal sends i — > —i, hence the requirement 
that the bosonic wavefunction is real implies that 

*L m (xi,X 2 , • • • ) = J] e -2^g(x i ,x 3 .) ^ ferm ( XljX2! . . . ) . 

(47) 

Thus, the time-reversal invariance of the bosonic wave func- 
tion is a highly nontrivial condition for the fermionic wave- 
function. 

For this reason, it is convenient to define a modified time- 
reversal transformation, which corresponds to naive time re- 
versal for the lattice fermions: 

2ferm : rfx -* dx! a -* -a; e -> e; i — > -i. (48) 

In first-quantized language 7f crm corresponds to complex con- 
jugation of the wavefunction for the fermionized vortices, 
^fcrm — * ^fcrm- It is important to emphasize that this trans- 
formation is not a symmetry of the fermionic Hamiltonian, 
since under 7f erm the field a changes sign while the Chern- 
Simons field A from Eq. g3 remains unchanged. 

However, as we shall argue further below, since the vor- 
tices interact logarithmically, their density fluctuations are 



greatly suppressed, and it is plausible that the phase factors in 
Eqs. J45i and ^47^ might be essentially the same for all vortex 
configurations that carry substantial weight. For an extreme 
example, if the vortices form a perfect charge-ordered state, 
then the Chern-Simons phase factor is constant. We conjec- 
ture below that this is also the case at the critical points sep- 
arating the spin-ordered states from the spin paramagnet. In 
any event, in situations where the Chern-Simons phase fac- 
tors are roughly constant for all important configurations, the 
physical and modified time-reversal transformations become 
essentially identical. It is then legitimate to require that the 
theory be invariant under 7f crm . This will be useful when we 
describe a conjectured low-energy continuum fermionic the- 
ory for criticality in the vortex system, since 7f erm acts in a 
simple, local way on the continuum fermion fields. 

B. Naive Continuum Theory 

To arrive at a low-energy continuum theory, we first con- 
sider a "flux-smeared" mean-field state with A xx > = 0. We 
will also ignore fluctuations in o xx ' for the moment, taking 
flxx' = 0. We are then left with half-filled fermions hopping 
on the honeycomb lattice with no fluxes. 

Diagonalizing the hopping Hamiltonian in momentum 
space using the two-site unit cell shown in Fig. [3] one finds 
that there are two Dirac points at momenta ±Q, where Q = 
(47r/3, 0). (Note that these are the same wavevectors found 
for the low-energy spin-1 excitations in the continuum analy- 
sis of the original spin model in Sec. Ill BP Focusing only on 
low-energy excitations in the vicinity of these momenta, the 
fermion operators can be expanded around the Dirac points. 
Denoting the lattice fermion field on the two sublattices as 
d xa with a = 1, 2 the sublattice label, we write 

to « ^flae lQ ' x + ialpTpLpe-^*, (49) 

where x continues to label the real-space position of the hon- 
eycomb lattice site and a y is a Pauli matrix. The fields ipR/L 
are two-component spinors that vary slowly on the lattice 
scale. Using this expansion we obtain for the free-fermion 
part of the continuum Hamiltonian 

Hf ] = [dx^v (p x a x + Vv ° v )^a , (50) 

where p = — i V is the momentum operator, an implicit sum- 
mation over the flavor index a £ [R, L) is understood, and we 
have suppressed the summation over the spinor indices a, (3. 
From now on, we absorb the nodal velocity v in the scaling of 
the coordinates. 

It is convenient to work in the Euclidean path integral de- 
scription. The free-fermion Lagrangian density is written in 
the form 

= ^a^d^a, (51) 

^R/l = tfn/L-f, (52) 

with implicit sums over the flavor index a £ [R, L] and space- 
time index /i £ [0, 1, 2], defined so that 90.1,2 = d T . x . y . The 
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2x2 Dirac matrices 7^ are given by 7 = a z , 7 = a y , 
7 2 = — <j x and satisfy the usual algebra {7^,7"} = 2<J M ". 
These matrices act within each two-component field so that 
(7 M V'a)a = lapipap- We will also find it useful to define 
Pauli matrices t x,v,z that act on the R/L flavor indices, i.e., 

Including the gauge-field fluctuations and the short-range 
fermion interactions, we arrive at a Lagrangian of the form 

1 i 
+ -^(e^xd^ax) 2 + —e^xA^Ax + Af ■ (53) 

Here e^x is the anti-symmetric tensor, and for simplic- 
ity we wrote a space-time isotropic form for the Maxwell 
action of the gauge field a M . We also used the standard 
(2 + l)-dimensional form for the Chern-Simons action that 
ensures the attachment of 2ir flux to the fermions, restoring 
the bosonic exchange statistics of the vortices. In the absence 
of the four-fermion terms C^j, there is a global SU(2) flavor 
symmetry, whose action on the fermion fields is generated by 
the r k Pauli matrices. 

The four-fermion interaction terms can be written as 

U S = A Jl + \ z Jl + A X (JJ + Jy) + \[3al^a)\ (54) 

where we have defined a flavor SU (2) vector of fermion bi- 
linears, 

J = ip a T ab ip b = ipla z T ab ip b , (55) 
and a "mass term," 

^0 = V'aV'a = ^l<J Z 1p a . (56) 

The four-fermion terms arise from vortex density-density in- 
teractions and other short-range interaction processes. For ex- 
ample, the Aq term roughly represents an overall vortex repul- 
sion, while the X z term represents the difference in repulsion 
between vortices on the same sublattice and opposite sublat- 
tices of the honeycomb. £4 / contains all independent four- 
fermion terms that can arise from the microscopic short-range 
fermion interactions (including possible short-range pieces 
mediated by the gauge fields). 

The above constitutes the naive continuum limit obtained 
by inserting the slow-field expansion Eq. ( I49> into the mi- 
croscopic Hamiltonian and assuming small fluctuations in the 
gauge field a and the Chern-Simons field A. We now dis- 
cuss the symmetries of the continuum formulation. Table II Vl 
shows the transformation properties of the continuum fermion 
fields deduced from the lattice fermion transformations in Ta- 
ble [ni| Also shown are the symmetry transformation prop- 
erties of the fermion bilinears J and Jo. Missing from the 
table is the original time-reversal invariance, which we do not 
know how to realize in the continuum. In its place we have 
included the modified time reversal 7f crm , defined in Eq. J48> . 
which corresponds to the naive time-reversal transformation 
for fermions hopping on the honeycomb lattice. Remarkably, 
the transformation properties of the flavor SU{2) vector J of 



TABLE IV: Transformation properties of the continuum fermion 
fields tpn/L an d the fermion bilinears J, Jo under the microscopic 
symmetries. Here J± — J x ± ij y . The spatial coordinates, which 
should transform appropriately under the lattice symmetries, have 
been suppressed for conciseness. The last column does not corre- 
spond to the physical time-reversal symmetry, which we cannot real- 
ize in a local manner in the continuum theory. Rather, it corresponds 
to the naive time reversal in Eq. 1481 for fermions hopping on the 
honeycomb lattice. 
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fermionic bilinears are identical to the transformation proper- 
ties of the bosonic flavor SU(2) vector defined in Sec. Ill Bl 
K — QaTab&b- Surprisingly, the fermionic bilinear J trans- 
forms under the modified time-reversal transformation 7f crm 
in precisely the way that the bosonic bilinear K transforms 
under the physical time reversal T. 

One can verify that the first four symmetries in Table IIVI 
(translation, rotation, modified reflection, and particle-hole) 
preclude all quartic fermion terms from the Lagrangian except 
the four terms exhibited in Eq. ( 154^ Moreover, these symme- 
tries prohibit all fermionic bilinears in the Lagrangian except 
a mass term of the form MJq = Mipip. This merits some 
discussion. In situations where it is legitimate to replace the 
time-reversal transformation by the modified symmetry Tf erm , 
we can use this modified symmetry to preclude such a mass 
term. However, generally this might not be possible. Indeed, 
as we shall see in the next subsection, a large mass term of this 
form when added to the Lagrangian places the system in the 
spin paramagnetic phase, provided the mass M has a specific 
sign relative to the Chern-Simons term. Since we know that 
the spin paramagnet does not break time-reversal symmetry 
T, it is clear that at least in this phase it is not legitimate to 
replace T with 7f crm . The reasons for this will become clear 
in the next subsection where we describe how the spin param- 
agnet and the spin ordered phases can be correctly described 
using the fermionized-vortex formulation. The discussion of 
whether or not it is legitimate to replace T with 7f crm at the 
critical points separating these phases will be deferred until 
Sec.UVDl 



C. Phases in the Fermionic Representation 

We now discuss how to recover the phases of the original 
spin model using the fermionic vortex fields. This extends our 
earlier discussion using the bosonic vortices in Sec. lIIIBl The 
paramagnetic phase of the original spin model is the bosonic 
vortex superfluid (Higgs) phase. In terms of the fermionic vor- 
tex fields, the paramagnet is obtained as an integer quantum 
Hall state, which in the continuum description corresponds 
to the presence of a mass term Mipip. Both fermion fields ipR 
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and ipL then have the same mass M with the same sign. Let us 
see that we indeed recover the correct description of the orig- 
inal spin paramagnet. Integrating out the massive fermions 
induces a Chern-Simons term for the sum field (a + A)^, so 
the Lagrangian for the fields a M and becomes 

1 i 
£ a ,A = ■^{e^\d l ,a\) 2 + —e^vxA^dvAx (57) 

+ zsign(M) g + A + (5g) 

47T 

We now stipulate that the sign of the mass M is taken to can- 
cel the original Chern-Simons term for the field A. One can 
then verify that the spectrum corresponding to the above La- 
grangian is gapped. For example, upon integrating out A the 
gauge field a obtains a mass. This is as expected, since there 
are no gapless excitations in this phase. Gapped quasiparticles 
of the paramagnetic phase of the spin model are described as 
follows. Consider acting with the fermion field ^ R j h on the 
ground state. This is essentially the same as inserting a vortex 
in the vortex field 9. The added fermion couples to a + A, 
and by examining £ a ,A we find that it binds V x a = —2ir 
flux to the fermion. Thus, the fermion is turned into a local- 
ized bosonic excitation, which is the familiar screened vortex 
in the vortex field 6 (the sign of the flux is of course consistent 
with our minimal coupling convention V6* + a). In terms of 
the original spin model, this bosonic excitation carries spin 1 
and is precisely the magnon of the paramagnet. 

It is worth remarking about the role of time-reversal sym- 
metry and the modified time-reversal transformation 7f erm of 
Eq. J48I in the spin paramagnetic phase. Since the spin para- 
magnet corresponds to an integer quantum Hall state for the 
fermions, it is clear that 7f erm will not be respected in this 
phase. This is consistent with Table IIVI which shows that 
the mass term ipip is odd under 7f 0lm . On the other hand, the 
phase factors in the fermionic integer quantum Hall wavefunc- 
tion will conspire to cancel the Chern-Simons phase factors in 
Eq. (147 \ leading to a wavefunction for the bosonic vortices 
which is real — consistent with physical time-reversal invari- 
ance. 

Consider now the magnetically ordered phases of the origi- 
nal spin model. These correspond to vortex insulators and are 
obtained in the fermionic theory as a result of spontaneously 
generating a fermion mass of the form 

(J) = (?H>) £ . (59) 

In the presence of such mass terms we have two massive 
Dirac fermion fields with opposite-sign masses, and integrat- 
ing out the fermions produces only a generic Maxwell term 
for the sum field (a + A)^. The gapless photon mode of 
the gauge field a then corresponds to the spin wave of the 
magnetically-ordered phase. Acting with a fermion creation 
operator also binds 2ir flux of the Chern-Simons field A and 
turns the fermion back into the original bosonic vortex. Be- 
cause of the gapless gauge field a such an isolated vortex costs 
a logarithmically-large energy as expected in the spin-ordered 
phase. 

Spontaneous mass generation is driven by interaction terms 
in the Hamiltonian such as the X z and Aj^ interactions. The 



(a) 

Vortex CDW 




FIG. 4: Vortex insulators that correspond to the spin-ordered states 
of Fig. Q (a) Staggered vortex "charge density wave" (CDW) state, 
corresponding to the coplanar v3 x state in the original spin vari- 
ables. Vortices preferentially occupy one of the two sublattices. (b) 
"Plaquette" vortex valence bond solid (VBS), corresponding to the 
dice-lattice collinear antiferromagnetic spin order. Vortex density is 
spread over the bold hexagons, (c) Vortex VBS, corresponding to the 
hexagonal collinear antiferromagnetic spin order. Vortices preferen- 
tially occupy the lattice of bold links. 



details of the magnetic order are determined by the specific 
mass term that is generated. For example, the mass term m z J z 
in the Lagrangian corresponds to the vortex "charge density 
wave" (CDW) state wherein vortices preferentially occupy 
one sublattice of the honeycomb lattice as shown in Fig. |4t. 
Indeed, from Table llVI this mass term is odd under the 7r/3 ro- 
tation and particle-hole symmetries, and can be identified with 
a staggered chemical potential for the lattice fermions that se- 
lects one of the charge-ordered states over the other. Therefore 
such spontaneously-generated mass gives rise to the transla- 
tion symmetry breaking in the vortex system that produces 
the CDW state. In the original spin model, this corresponds to 
the coplanar spin state of Fig. 

On the other hand, the mass term m x J x + m y J y corre- 
sponds to the vortex valence bond solid (VBS) state. The 
specific pattern is not resolved with only the four-fermion in- 
teractions, but we expect that because of the underlying lattice 
there are higher-order terms in the action that pin the direction 
in the (J x , J y ) plane so that (J+) 3 is either +1 or — 1, The 
corresponding bond orders are deduced by interpreting the 
spontaneously-generated mass term \m\(e %a J+ + e~ la J-)/2 
as inducing a modulated vortex hopping amplitude t x x > = 
1 + \m\ cos[Q • (x + x') + a]/2. For e i3a = -1 the stronger 
bonds form hexagons and produce the vortex "plaquette" VBS 
shown in Fig. [4j% which in the original spin model corre- 
sponds to the dice-lattice collinear spin order of Fig.^. On 
the other hand, for e l3a = 1 the strong bonds form the lattice 
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shown in Fig.|4j;; the resulting vortex VBS corresponds to the 
hexagonal collinear spin state in Fig.[Q;. 

It is worth recalling from Sec. lIIBI that in the Landau theory 
of the spin model, the chirality order parameter (k) that devel- 
ops in the coplanar state and the bond energy wave (complex) 
order parameter (B) that is non-vanishing in the collinear spin 
states are both expressible in terms of the SU (2) vector K 
defined in Eq. (1221 : 

k~(K z ); B~(K x + iK y ). (60) 

Thus, the magnetically-ordered states are characterized by 

(K) = ($*r$) + 0, (61) 

with coplanar order corresponding to ordering K z and 
collinear order corresponding to an ordering in the (K x , K y ) 
plane. This is directly analogous to the ordering patterns of 
the SU(2) vector J = tpTi/j as discussed above. 

The fermionic formulation of the dual vortex theory thus 
allows us to correctly describe the phases of the original spin 
model. From the point of view of the bosonic vortex sys- 
tem at half-filling, this is rather nontrivial. For example, the 
fermionic formulation correctly captures the two low-energy 
spin-1 excitations with wave vectors ±Q in the spin param- 
agnet phase, even though these wavevectors are not appar- 
ent in terms of the bosonic vortices. Moreover, the magnet- 
ically ordered phases are accessed in a unified manner within 
the fermionic formulation, via a spontaneous mass generation 
driven by the vortex interactions which produces either the 
vortex CDW or one of the VBS states, thereby "releasing" 
the dual photon. While the vortex CDW is perhaps natural 
in the bosonic vortex theory given the strong repulsive vortex 
interactions and the specific charge ordering, the VBS phases 
are nontrivial for the vortex bosons at half-filling. Indeed, a 
common approach for analyzing such VBS states of bosons 
is to study their dual formulation, 36 which in the present con- 
text corresponds to our analysis of the original spin model. 
Our alternative route via the fermionization achieves this due 
to the fact that it in some sense combines the direct and dual 
perspectives. 

D. Criticality in the Fermionic Theory? 

Encouraged by these successes, we now embark on a study 
of the transitions between the phases using the fermionic for- 
mulation. As a guide we use the anticipation that the original 
spin model has continuous transitions. We have confidence 
for the presence of the 0(4) fixed point and its likely appear- 
ance as a multicritical point, and also for the continuous de- 
coupled U(l) x £7(1) transition between the paramagnet and 
the collinearly-ordered phase. Moreover, the paramagnet-to- 
coplanar ordering transition may also be continuous. We want 
to now see if we can access any of these critical points within 
the Dirac theory of fermionized vortices. 

At criticality, we expect the vortices to be massless. We also 
expect the original spins to be critical, which corresponds to 
the dual gauge field a being critical. To proceed we first per- 
form the following trick in the continuum action. The fermion 



current couples to the combination 5^ = 0^+ A^, and it is 
instructive to retain only this field in the path integral by per- 
forming the Gaussian integration over the field A. The result- 
ing naive continuum Lagrangian has the form 

£ = V^VC^ - ia^tpa + £4/ 

1 7T 

+ 7-^(V x a) 2 + i—r{X7 x a) ■ (V x V x a) . (62) 

2e z e 4 

Since a should be massless at criticality, it seems likely 
that the last term above (which can be viewed as a higher- 
derivative Chern-Simons term) is irrelevant compared with the 
Maxwell term at long wavelengths. Indeed, ignoring the ef- 
fects of this term on the field gives a correct description in 
the vortex insulator phases, where the fluctuations of are 
essentially the same as a M and represent gapless spin waves 
of the original spin model. As a putative description of criti- 
cality in the original spin model, we henceforth drop this term 
altogether. The theory is then equivalent to quantum electro- 
dynamics in (2+1) dimensions, 

Cqeds = ?„7"(fy. - *a M )Va + ^(Vxi) 2 + £ 4J .(63) 

We conjecture that this QED3 theory with two Dirac fermions 
contains a description of the critical properties of the original 
spin model. 

Words of caution are in order here. Dropping the higher- 
derivative Chern-Simons term changes the symmetry of the 
continuum Lagrangian, which is now invariant under the mod- 
ified time-reversal symmetry 7f crm , as can be seen from Ta- 
ble^] Apparently, neglecting the higher-derivative Chern- 
Simons term is tantamount to replacing the physical time- 
reversal symmetry transformation by the modified one. This 
seems reasonable on physical grounds since the vortices are 
strongly-interacting, and power-counting in the Lagrangian 
lends further mathematical support to the validity of this ap- 
proximation. Once we have replaced T by 7f erm , both the 
fermion mass term Mtptf) and the higher-derivative Chern- 
Simons term are symmetry precluded. This seems consistent 
since we are looking for a fixed-point theory with massless 
fields. 

It is worth emphasizing that the above QED3 Lagrangian 
is the proper continuum theory for fermions hopping on the 
honeycomb lattice coupled to a noncompact gauge field (no 
Chern-Simons field). Invariance under the modified time- 
reversal symmetry 7f crm follows provided the fermionic hop- 
ping Hamiltonian is real. Our proposal is that such a critical 
QED3 theory faithfully describes the continuous transitions 
of our bosonic vortices at half-filling. Again, the strong log- 
arithmic interaction between vortices greatly suppresses their 
density fluctuations, and one expects that as a result the vor- 
tex statistics might not be so important. The trick that elimi- 
nated the Chern-Simons gauge field A leaving only the higher- 
derivative Chern-Simons-like term for a can be viewed as a 
more formal argument for this. If we are allowed to drop such 
higher-derivative terms, then we can essentially eliminate the 
statistical Chern-Simons field by absorbing it into the already- 
present gauge field a, which is precisely the field that mediates 
the vortex repulsion. 
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Assuming that QED3 is appropriate for describing critical- 
ity in the spin model, it is not at first clear whether or not this 
requires a further fine-tuning of the mass term Mipip to zero in 
the continuum theory. However, once we assume that it is le- 
gitimate to replace T by 7f erm at criticality, then this fermion 
mass term is symmetry-precluded, and we are not allowed to 
add it to the critical Lagrangian as a perturbing field. 

In the next section we will analyze the QED3 Lagrangian 
per se, focusing on its potential to describe criticality in the 
spin model. 



V. ANALYSIS OF QED3 

The QED3 theory with two Dirac fermions as realized on 
the half-filled honeycomb lattice is a difficult problem with 
its own phase diagram. As we will argue below, it is likely 
that the lattice model generically ends up in a phase with a 
spontaneously-generated fermion mass. In this case, the con- 
tinuum Lagrangian Cqed3 with massless fields potentially 
applies only to critical points of the lattice model. 



A. Critical Theory 

In addition to the discrete symmetries tabulated in Table llVI 
the Lagrangian Cqed3 has a number of continuous global 
symmetries. The full Lagrangian is invariant under the (dual) 
global U(l) symmetry, ip a — > e la ip a - Associated with this 
symmetry is a conserved vortex current G M = tpj^ip, which 
satisfies d^G^ = 0. Here G° is simply the vortex density. 
In the absence of the four-fermion interaction terms, the La- 
grangian also enjoys a global flavor SU(2) symmetry, being 
invariant under 

il> a -> U ab ip b ; tp a -> ip b Ul a , (64) 

for arbitrary £17(2) rotation, W — U —1 . The associated 
SU(2) conserved currents are given by j M = ip^Tip, and 
satisfy d^j^ = 0. The remaining fermionic bilinears Jo = 
ifrtp and J = tpT'ip are not parts of any conserved current; Jo is 
a scalar while J rotate as a vector under the flavor SU{2). In 
the special case \ z — Aj_, the four-fermion interaction terms 
are also invariant under SU(2) flavor rotations, but more gen- 
erally are not. 

Due to the gauge interactions QED3 is a strongly- 
interacting field theory 37 . Specifically, expanding about the 
free-field limit with e 2 = and £4/ = 0, the continuum 
fermion fields have scaling dimension A^ = 1 so that the 
four-fermion interactions are irrelevant. However, gauge in- 
variance dictates that a has scaling dimension Aa = 1, im- 
plying that e 2 is relevant and grows in the infrared. 

To seek a controlled limit one is forced to generalize the 
model in some way. Perhaps the simplest approach is to in- 
troduce N copies of the fermion fields, t/j aa — > ipj aa with 
j = 1,...,N, each of which is minimally coupled to the 
same gauge field, and to then study the model in the large- 
ly limit, 37 ] 38 ! 39 ! 40 ] 41 ] 42 (Note that the theory with two fermion 



flavors corresponds to N = 1, and the flavor symmetry of the 
theory thus generalized is SU (2N).) Upon integrating out the 
fermions and expanding to quadratic order in the gauge field 
one obtains an effective gauge action of the form, 

Sa = Tglwf^ " W) ~ a ^- q) ~ a » {q) ' (65) 

with a small coupling constant g = 8/N. The gauge prop- 
agator proportional to g/\q\ mediates a screened interaction 
between the fermions which falls off as l/\r\ in real space 
and is much weaker than the bare logarithmic interaction. At 
infinite N the gauge fluctuations are completely suppressed, 
and, except for some subtleties that will not be important 
here^ 8 the model reduces to free Dirac fermions. It is then 
possibla 3 i 3 ^22i 4 iii 4 ii2 2 i to perform a controlled analysis pertur- 
bative in inverse powers of N. Specifically, one can compute 
the scaling dimension of various perturbations, such as the 
quartic fermion terms, order by order in 1/N. 

To obtain the leading 1/N corrections it suffices to re- 
tain only the original two fermion fields ipR/L that appear 
in Cqed3 and to replace the Maxwell term by the singular 
gauge interaction from Eq. J65I . One can then perform a sim- 
ple Wilsonian renormalization group analysis perturbative in 
the single coupling constant g. After integrating out a shell of 
modes between a cutoff A and A/6 with b > 1, the fermion 
fields can be rescaled to keep the tpj^daip term unchanged. 
Gauge invariance then automatically ensures that ^pj^a^tp is 
also unchanged. Due to the singular momentum dependence 
in Eq. (I65l >. 1/g cannot pick up any diagrammatic contribu- 
tions from the high-momenta mode integration. With Ag = 1 
assured by gauge invariance, rescaling will not modify g, and 
the theory describes a fixed line parameterized by the coupling 
9- 

With this simple Wilsonian renormalization group in hand, 
one can easily compute the scaling dimensions of the quartic 
fermion operators perturbatively in g. We find that to leading 
order in g the scaling dimension of the X' term is unmodified, 
withA^ =4 + 0(.g 2 ). 

The other three quartic fermion terms mix already at first- 
order in g. Of the three renormalization group eigenoperators, 
we find two that are singlets under the flavor SU (2), 

Qo = J 2 - Jo ; Q'o — J 2 + 3Jo , (66) 
and one that transforms as spin 2 under the flavor SU (2), 

Q 2 = J 2 - 3J 2 . (67) 
To first order in g the respective scaling dimensions are 

A Qo =4-4 ff /7r 2 ; Aq' o = Aq 2 =4 + Ag/ (3ir 2 ). (68) 

The above discussion is based on a particular generaliza- 
tion of the terms in £4/ to the SU (2 N) -symmetric theory, 
where quartic terms contain only two fermion flavors. An- 
other natural generalization proceeds by classifying all four- 
fermion terms in the SU (2A)-symmetric theory according to 
the irreducible representation of flavor and Lorentz symme- 
try under which they transform. It is possible to establish a 
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natural correspondence between four of these multiplets and 
the multiplets in the SU (2) theory to which the terms in £4 1 
belong. One can then calculate the scaling dimensions of the 
resulting terms order-by-order in 1/N. For the terms corre- 
sponding to Qo, Q' and Q2, we reproduce the same scaling 
dimensions above. On the other hand, the analog of the A 
term has dimension 6 + 0(1/N), as can be seen by a calcu- 
lation of its autocorrelation function at N = 00. We remark 
that both generalizations of this term strongly suggest it is an 
irrelevant perturbation. 

At sufficiently large N all of the quartic terms have scal- 
ing dimensions greater than D = 2 + 1, and are thus irrel- 
evant. The Lagrangian Cqed3 is then critical and describes 
a conformally-invariant, strongly-interacting fixed point. Our 
hope is that this fixed point (at N = 1) corresponds to one of 
the three critical points of the original spin model discussed in 
Sec. Ill Dl We will try to identify which one in the next subsec- 
tion. In Sec. IV CI we will return to the important issue of the 
stability of this fixed point in the physically-relevant N = 1 
case. 



B. Multicritical Point in the Spin Model? 

In the absence of the quartic terms, the critical theory de- 
scribed by Cqed3 has a global SU(2) symmetry shown ex- 
plicitly in Eq. (1641 . Since the three fermionic bilinears com- 
prising J rotate as a vector under this SU(2) symmetry, at 
the critical point each component must have the same scal- 
ing dimension. As noted in the previous section, the vector J 
and the vector of bosonic bilinears K have identical symme- 
try properties under all of the microscopic lattice and internal 
symmetries, and can thus be equated at criticality. So we then 
ask: At which of the three putative critical points in the origi- 
nal spin model does K transform as an SU(2) vector? 

Recall that the transition between the paramagnet and the 
collinearly-ordered spin state is described by two independent 
£7(1) models, one for each of the two complex fields §r/l- 
Since K + = §* r $l, it is evident that K z = \§r\ 2 - |$l| 2 
most definitely does not have the same scaling dimension as 
K x and K y at this transition. Moreover, at the transition from 
the paramagnet into the coplanar spin state with long-range 
chiral order, one expects the chirality k ~ K z correlator to 
decay substantially more slowly than the bond density wave 
operator B ~ K + , which only orders with collinear spin or- 
der. This leaves only the multicritical point in the spin model 
as a candidate fixed point described by QED3 with the full 
SU{2). 

Since the multicritical point is the special point where the 
paramagnet merges with both the coplanar and the collinear 
spin-ordered states, one expects both K z and K± to have 
slowly-decaying correlators. In fact, as the Landau-Ginzburg- 
Wilson analysis demonstrates, the multicritical point has an 
emergent global 0(4) symmetry. Moreover, the three compo- 
nents of K together with six other bosonic bilinears transform 
as a symmetric and traceless second-rank tensor representa- 
tion of 0(4). Indeed, the vector K transforms as a vector 
under an SU(2) subgroup of the 0(4). We are thus led to the 



rather bold conjecture: The critical theory described by QED3 
with two Dime spinors is identical to the 0(4) critical point 
o/$ 4 scalar field theory. 

As detailed in Sec. [II] at the 0(4) critical point of scalar 
field theory, the operators can be arranged into multiplets 
which transform under irreducible representations of the 0(4) 
symmetry group. For example, the real and imaginary parts of 
Qr/Li denoted \i with i — 1, . . . , 4 in Eq. (1241 . transform as a 
vector under 0(4), while the ten independent bilinears formed 
from XiXj decompose into an 0(4) scalar, Kq — Yli=i Xi> 
and the nine components of a symmetric, traceless matrix that 
transforms as a second-rank tensor under 0(4). Moreover, at 
the 0(4) multicritical point of the spin model there are pre- 
cisely two relevant operators that must be tuned to zero: the 
quadratic mass term 

Qo = K = \<$> R \ 2 + |$ L | 2 , (69) 

which is an 0(4) scalar, and the quartic term 

Q 2 = K 2 - 3K 2 , (70) 

which transforms as spin-2 under 0(4) and breaks the global 
0(4) symmetry down to £7(1) x £7(1). When the coupling of 
Q2 is positive in the Hamiltonian, it favors the coplanar spin- 
ordered state {(K z ) ^ 0) over the collinear states ((K XiV ) ^ 
0), and vice versa when it is negative. 

In order to back up our bold conjecture that the critical 
QED3 theory describes this 0(4) multicritical point, it is 
clearly necessary, at the very least, to find operators in the 
QED3 theory that can be associated with each of the 0(4) 
multiplets mentioned above. That is, one must find opera- 
tors in the QED3 theory which under the microscopic sym- 
metries transform identically to their Landau theory counter- 
parts. Moreover, one should identify the two relevant oper- 
ators in the QED3 theory that transform identically to Qo,2 
above. Completing the latter task requires revisiting the issue 
of stability of the QED3 fixed point. 

C. Stability of the QED3 Fixed Point 

To assess the stability of the QED3 fixed point, we need to 
determine which of the quartic interactions are relevant per- 
turbations. Equivalently, we need to see which of the quartic 
interactions have scaling dimensions smaller than D — 3. In 
general this is a very difficult task. The best we can do is to 
examine the trends in the 1 /N expansion. To leading order 
in 1/N only the SU(2) scalar Qo in Eq. d66> has a scaling 
dimension which is reduced below 4, and if we naively put 
N = 1 we find Aq = 0.76 < 3. Let's assume that Qo is 
indeed relevant at the N = 1 critical QED3 fixed point. Can 
we then identify Qo with one of the two relevant perturba- 
tions at the 0(4) Landau critical point? Since Q is invariant 
under all of the symmetries, it has the same symmetry as the 
Landau mass term, Q = K . Moreover, if we add Q to 
the fixed-point QED3 Lagrangian with a positive coupling the 
paramagnetic spin state with (Jo) 7^ is favored, whereas a 
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negative coupling drives magnetic spin order, (J) 7^ 0. Like- 
wise, a positive mass term Qq when added to the 0(4) sym- 
metric fixed point of the Landau action drives one into the 
paramagnet, while a negative mass leads to magnetic order. 
Apparently, it is then entirely consistent to identify the mass 
term Q in the Landau theory with the QED3 perturbation Q . 

Similarly, it is possible to identify the second relevant per- 
turbation at the 0(4) Landau fixed point, Q 2 = K 2 — 32C^, 
with the quartic perturbation Q2 = J 2 — "iJ 2 in QED3. We 
have already established that the two vectors K and J trans- 
form identically under all of the microscopic symmetries (see 
tables HI and lIV> . Moreover, adding the two operators to their 
respective Lagrangians breaks the energetic degeneracy be- 
tween the coplanar and collinear spin-ordered phases. If this 
identification is indeed correct, it would imply that the quar- 
tic perturbation Q 2 is (weakly) relevant when added to the 
critical QED3 Lagrangian. Based on our leading-order 1/N 
calculation this is surprising since we found that the lead- 
ing correction increased the scaling dimension of this oper- 
ator: Aq 2 = 4 + 4g/(3n 2 ) + 0(g 2 ). But it is conceivable 
that higher-order corrections are negative and dominate when 
N = 1. 

At this stage our precise operator correspondence between 
0(4) scalar field theory and QED3 has been limited to op- 
erators which are invariant under the global 17(1) spin sym- 
metries. In order to complete the task, then, we must iden- 
tify operators within QED3 that are symmetry-equivalent to 
the spin fields, &r/l, as well as the anomalous bilinears, 
I = &a[T TV }ab&b an d I ■ Since these operators carry non- 
zero spin, they correspond to operators in QED3 which add 
(dual) gauge flux (A x a) in discrete units of 2tt. We study 
the properties of such "monopole" operators in the next sub- 
section. 



D. Monopole Operators in QED3 

It is important to note that since we are dealing with a non- 
compact gauge theory, we do not have monopoles as dynam- 
ical degrees of freedom in the system. We can still, how- 
ever, consider states with monopoles inserted by hand. Our 
treatment is motivated by the work of Borokhov et a/. 43 , 
who discussed such topological disorder operators in three- 
dimensional QED. 

Let us first consider the monopole operators in our dual 
QED3 formulation that correspond to &r/l an d &*r/l- Since 
these fields carry spin ±1, the corresponding monopole op- 
erators add ±2ir gauge flux. In what follows, we will treat 
this gauge flux as a classical, static background to which the 
fermions respond. Moreover, we will assume that the flux is 
spread out over a large area compared to the lattice unit cell 
since these configurations have low energy. In the presence of 
such monopoles, the Hamiltonian density in Eq. (I50i general- 
izes to 

nf q = V4 [(Px - a£K + ( Py - a«)a*]Va , (71) 
where as before a = R/L denotes the fermion flavor and is 



implicitly summed. In Eq. M\\ . a q x denotes the vector po- 
tential giving rise to 2irq flux, where q E Z is the monopole 
"charge". We will define monopole creation and annihilation 
operators by specifying their action on the zero-flux state, and 
therefore focus on the q = ±1 sectors. 

In the presence of a charge q = ±1 monopole, the above 
Hamiltonian admits two zero-energy states, one for each 
fermion flavor. 44 Moreover, one of these states must be oc- 
cupied in order to obtain a physical (i.e., gauge-invariant) 
q = ±1 ground state. To see this explicitly, first note that 
since the honeycomb lattice is bipartite, sending d x — > — e? x 
on one of the two sublattices changes the sign of the fermion 
hopping part of the lattice Hamiltonian, Eq. (14-1 1 . Hence there 
is a one-to-one correspondence between states with energy E 
and — E. Since the fermions are at half-filling, the lowest- 
energy q = ±1 states will therefore have all of the E < 
modes and one of the two zero modes occupied. 

Thus, we see that our dual QED3 by its own dynamics has 
two species of monopole insertions carrying R/L momentum 
(see below). Such topological disorder operators are expected 
to have nontrivial scaling dimensions in the interacting con- 
formal field theory4£ We can also ask how these monopoles 
transform under the various microscopic symmetries. For the 
lattice symmetries, we can guess the answer from the trans- 
formation properties of the i/jr/l fermion fields in Table Hvl 
by simply retaining the t matrix structure acting on the R/L 
flavor indices. However, the transformations of the fermion 
fields have U(l) ambiguity, while the monopole insertions are 
gauge-invariant objects and therefore have definite transfor- 
mation properties. To fix this ambiguity we need to examine 
the monopole states more carefully, which is outlined below. 

The analysis is performed in a fixed gauge and starts by 
considering the zero modes and how they transform under the 
microscopic symmetries. We obtain the zero-energy states of 
Eq. (17 1 i in the presence of a charge q = ±1 monopole fol- 
lowing the approach of Jackiw in Ref. 44. In the Coulomb 
gauge, the vector potential may be written in terms of a scalar 
function a q as af = eijdja q . Assuming an azimuthally- 
symmetric flux distribution centered around the origin, we can 
choose a q ~ — gln|x| as x — > 00 to achieve a total gauge 
flux of 2nq. Focusing only on the unit-strength monopoles, 
we make the replacement 



V'a(x) -> <j>aq{x)faq 



(72) 



in Eq. M\\ . where for each a = R/L and q — ±1, Q9 (x) is 
the corresponding two-component, zero-energy wavefunction 
and the operator f aq annihilates this state. It is straightforward 
to show 44 that with the above gauge choice, the only suitable 
wavefunctions are 



|x| (0 
WO 
x U 



(73) 
(74) 



These states are quasi-localized since their normalization in- 
tegrals diverge logarithmically with system size, but the only 
information that we actually need is the spinor structure of 
these wavefunctions. 
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TABLE V: Transformation properties of the zero-mode operators 
fn/L.q in the charge q — ±1 monopole sectors. Note that these 
transformations were obtained by employing the Coulomb gauge. 







Rtt/3 




C 


Tferm 




e tQdl fn, q 






J L.-q 


qh-q 


/i,9 —> 






fL,q 


J R,-q 


-qfR-q 



The symmetry transformation properties of the zero-mode 
operators f aq can be deduced from Eq. J72I together with the 
transformation properties of the continuum fields ip a in Table 
II VI Under translations for instance, ip R — ► e lC *' Sr tp R , which 
implies that f Rq — ► e lC *' Sr f Rq . Since the flux changes sign 
under particle-hole and fermionic time reversal, one must ad- 
ditionally transform q — > — q under these symmetries. Thus, 
particle-hole transforms ipR — > cr x [np L ] t and f R , q — > f* L , 
while time reversal sends ip R — * ia y ipL and f Rq — ► q]h,-q- 
The transformation properties of fn/L,q are summarized in 
Table|Vl 

With the zero modes in hand, we are now in position to 
introduce the monopole operators that correspond to §r/l 
and &* R / L - First, we define operators M R , L that insert +2ir 
flux and add momentum ±Q by filling the resulting R/L zero 
mode. Their Hermitian conjugates M r /l insert — 2ir flux and 
add momentum TQ by filling the L/R zero mode. In partic- 
ular, M R , L and M R / L act on the zero-flux ground state |0) 
according to 

M R/L \0) ee e ia ^f R/Lj+1 \DS,+l) (75) 
M R/L \0) ee e l ^^fl /R ^\DS,-l), (76) 

where \DS, ±1) denotes the filled negative-energy Dirac sea 
in the presence of ±2ir gauge flux, with both zero modes 
vacant. The right-hand sides of the above equations are 
gauge-invariant quantum states and can be chosen with ar- 
bitrary phases a R /L and (3 r /l (indeed, for any orthogonal 
|0),|X),|Y), we can define A ee \X){0\ + \0)(Y\ such 
that A\0) = \X) and A+jO) = \Y)). We emphasize that 
we have defined the monopole creation/annihilation opera- 
tors by constructing the corresponding physical states with 
one monopole/antimonopole, and we have the full freedom 
to specify the phases a R /L and f3 R /L to our convenience. 

The transformation properties of M r /l can be inferred 
from Eqs. (175 \ and M6\ . provided we know how the states 
|0) and \DS, ±1) transform. While |0) transforms simply, de- 
ducing the symmetry properties of the negative-energy Dirac 
sea in the presence of ±2tt flux is challenging since this would 
require the knowledge of all wavefunctions and not only the 
zero-energy modes. Nevertheless, even without such knowl- 
edge, by employing relations among the symmetries we can 
infer the needed transformation properties of \DS, q) so as 
to be able to fix the transformation properties of M R /^ al- 
most completely by appropriate choices of a R /^ and /3 r /l- 
As this is somewhat involved, we relegate the details to Ap- 
pendix[E] and simply recall the results here. We find that the 



filled negative-energy Dirac sea in the presence of a charge 
q = ±1 monopole transforms as follows: 

T 5r \DS,q) = \DS,q) 

R*/a\DS,q) = Cn/ 3 e-' l ^/ 3 \DS,q) 

n x \DS,q) = C*\DS,q) (77) 

C\DS,q) = fl_ q fl_ q \DS,-q) 

T [crm \DS,q) = \DS,-q). 

Here Ctt/3, Cx = ±1 are possible overall signs which have not 
been determined in our analysis. 

With the help of Table |VJ we now know the transformation 
properties of the gauge-invariant states in Eqs. d75l > and i76i . 
and it is easy to check that we can choose the phases a R /L and 
j3 R /L such that the monopole operators transform as follows, 



Tsr 


M R/L 


—>e ^ M R / L 


Rtt/3 


M R/L 


^M L/R 


Tlx 


M R/L 


^M R/L 


C 


M R/L 




^ferm 


M R/L 





Here the sign s = ±1 is in principle fixed, but can not be de- 
termined from our analysis. Comparing with TableQ]we see 
that we have constructed monopole insertion operators M r /l 
in the dual QED3 that transform identically to §r/l m the 
spin system under the lattice symmetries and particle-hole. 
Additionally, M R j L transform under 7f clm in the same way 
that Q r /l transform under the physical time reversal, up to 
a possible overall minus sign. We are thus content with the 
agreement obtained here between the discrete symmetries in 
the two systems (a more precise analysis would of course be 
welcome). 

We turn next to double-strength, charge q = ±2 
monopoles. Although the operators M R j L defined above are 
insufficient to create all such monopoles (see below), we can 
nevertheless deduce from M r /l the transformation properties 
of the most relevant double-strength monopole operators us- 
ing the following argument. For example, one can use M R ^ L 
to create two isolated q = +1 monopoles separated by a large 
distance. If we fuse these monopoles together, the quantum 
numbers of the resulting double-strength monopole will be 
the sum of those for the two unit-strength monopoles since 
the quantum numbers are conserved. Such double-strength 
monopole operators must therefore transform identically to 
the "bilinears" M R , M\, and M r Ml, where we treat M R / L 
as bosonic fields. These bilinears can be grouped into an 
5(7(2) vector, 

D = M a [T T y] ab M b , (79) 

which transforms identically to the SU (2) vector J defined in 
Eq. ( I26l l. Thus, we have established the correspondence be- 
tween the anomalous spin-field bilinears and double-strength 
monopole operators. 
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As remarked above, the preceding argument pertains only 
to the most relevant double-strength monopoles. In the pres- 
ence of 47r flux there are four zero modes (two for each fla- 
vor), and two of these need to be occupied to obtain a gauge- 
invariant state. Hence, one has to consider six double-strength 
monopole operators since there are six ways of filling the zero 
modes. However, three of these states are not locally gauge- 
neutral and have contributions to the gauge charge density 
with dipolar angular dependence. Such states have higher en- 
ergy due to associated electric fields than the remaining three 
states which have rotationally-invariant locally neutral charge 
distributions. Operators that create such dipolar states are ex- 
pected to have higher scaling-dimensions and are neglected, 
leaving only the three double-strength monopole operators 
considered above. 

To summarize our discussion, the monopole insertions 
Mr/l in QED3 correspond to the Qr/l fields in the spin 
model. We predict that the scaling dimension of such 
monopole operators in QED3 with N = 1 and full SU (2) 
flavor symmetry is given by the scaling dimension of the or- 
dering field $ at the 0(4) critical point of $ 4 scalar field the- 
ory. Furthermore, we identified double-strength monopoles 
D, in QED3 that correspond to the spin-field bilinears 
J, P . Recall also that the fermionic bilinears J correspond 
to the spin bilinears K. At the 0(4) fixed point, K, I, and V 
all have the same scaling dimension — these are the quadratic 
anisotropy fields of the $ 4 theory. Thus, we are led to a rather 
unusual prediction: the fermionic bilinears and the double- 
strength monopole insertions have the same scaling dimen- 
sion in our QED3 theory. We emphasize that this is striking 
from the perspective of the large- N limit, where the monopole 
operators have large scaling dimension, gauge fluctuations be- 
ing strongly suppressed, while the fermionic bilinears have 
scaling dimension 2 — 0(1/N). It would be quite remarkable 
if the scaling dimensions of these operators which have no ob- 
vious relation indeed merge in the physically-relevant N = 1 
case. 



E. Dual Picture for the Honeycomb QED3 

We conclude this section with a separate argument for the 
applicability of the critical QED3 theory to the description of 
the continuous phase transitions in the original spin model. 
We intend to further address the question of whether it is 
appropriate to drop the higher-derivative Chern-Simons-like 
term in the fermionized-vortex Lagrangian Eq. (1621 as irrele- 
vant at criticality and at the same time prohibit the mass term 
Mtpip. The setting below, which in many respects is a re- 
versal of the preceding line of attack (but not a circular ar- 
gument), gives us good control on the effects of the Chern- 
Simons gauge field in a situation where the particles in the 
theory already have strong interactions mediated by a photon 
field. Our difficulties in Sec. llVDl which we distilled to the 
question regarding the realization of the physical spin time- 
reversal symmetry in the continuum dual action, stem from 
the inability to treat the statistical interactions accurately in 
this context. Specifically, the local symmetries of the formal 



Lagrangian Eq. i62\ allow a generic mass term Mtpip which 
is of course a very relevant perturbation in the Dirac field the- 
ory. On the other hand, all such 7f olm -violating terms have to 
conspire to recover the original spin time-reversal invariance 
which appears lost in the continuum theory. 

We therefore look for a similar situation where the statis- 
tical interactions can be treated accurately. In the present 
context, it is natural to consider the honeycomb lattice QED3 
problem per se, and seek a "dual" description of this fermionic 
theory in a manner familiar from the treatment of "vortices" 
in fractional quantum Hall systems. 45 46 - 47 This is achieved by 
viewing the lattice fermions as bosons carrying fictitious 2ir 
flux and performing a duality transformation on these Chern- 
Simons bosons. The Chern-Simons bosons are governed by 
a Hamiltonian similar to Eq. d34i . but with additional statisti- 
cal interactions mediated by the corresponding Chern-Simons 
gauge field. The duality transformation essentially reverses 
the steps in Sec. |ni| but in the presence of the Chern-Simons 
gauge field coupled to the currents conjugate to the 9 fields. 
In the absence of the Chern-Simons gauge field, we would of 
course recover the original XY spin model on the frustrated 
triangular lattice. For example, in terms of the integer-valued 
currents j conjugate to the <fi fields, we would obtain a Eu- 
clidean action of the form 

SxyU]^^-^"'^, (80) 

where the first term denotes schematically generic short-range 
repulsion between the currents and the second term encodes 
the frustration for the ^-boson hopping. Including the Chern- 
Simons gauge field, the corresponding integration can be per- 
formed essentially exactly and leads to a new contribution of 
the form 

6S\j] = iKj2i-(y xj) (81) 

with some nonuniversal numerical coefficient k. Thus, the 
dual of the honeycomb lattice QED3 model is the stacked 
triangular lattice antiferromagnetic XY spin model with ad- 
ditional current-current interactions that break spin time- 
reversal and a; -reflection symmetry. 

The advantage here compared with our fermionized vortex 
system of Sec. lIVl is that the new interactions SS [j] are local 
in terms of the spin fields of the new spin model. We also ob- 
serve that the new model still respects the translation (Ts r ), ro- 
tation (R n / 3 ), modified reflection {1Z X ), and particle-hole (C) 
symmetries defined in Sec. Ill Bl Furthermore, the expected 
charge-ordered and integer quantum Hall fermionic states of 
the honeycomb QED3 model are recovered as, correspond- 
ingly, the spin-ordered and paramagnetic phases in the dual 
description. Performing the continuum analysis of the new 
spin model, we conclude that the above symmetries preclude 
all relevant terms other than the ones already exhibited in the 
action Eq. illl . In particular, the additional terms that violate 
spin time-reversal symmetry necessarily have more deriva- 
tives and cannot generate new relevant terms. We therefore 
conclude that these terms are irrelevant in the <f> 4 theories that 
govern criticalities in the spin model, and that the spin time- 
reversal symmetry emerges at the corresponding fixed points. 
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Thus, it appears that we have identified the continuum N = 1 
massless QED3 with the criticality in our original classical 
spin system discussed in detail in Sec.|ll] A possible flaw in 
this identification is that the new interactions SS[j] are nu- 
merically very large, and therefore put the spin system corre- 
sponding to the lattice QED3 into some very different regime 
of the phase diagram. We cannot rule this out with the present 
approach, but our predictions for the QED3 criticality can be 
in principle tested in direct lattice simulations. 



VI. CONCLUSION 



criticality of bosons at finite density in situations where the 
bosons are strongly-interacting, 20 - 50 leading to a diminished 
role of exchange statistics. 
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In this paper we have introduced a new dual vortex ap- 
proach to easy-plane frustrated quantum spin systems which 
is well-suited for the study of novel critical points and criti- 
cal phases. Specifically, we demonstrated how one can for- 
mulate a low-energy dual description of frustrated systems 
by fermionizing the vortices, which are at finite-density due 
to frustration, via Chern-Simons flux attachment. A detailed 
analysis was carried out for the easy-plane integer-spin tri- 
angular antiferromagnet, guided by the known results from a 
direct Landau analysis of the spin model. The fermionized- 
vortex approach led naturally in this case to a low-energy 
Dirac theory for two species of fermions coupled to a Chern- 
Simons field and a U(l) gauge field that mediates logarith- 
mic vortex interactions. We demonstrated how the ordered 
phases of the spin model are captured in the continuum dual 
theory via spontaneous fermion mass generation, and conjec- 
tured that the critical QED3 theory obtained by essentially ig- 
noring the Chern-Simons field describes the 0(4) multicrit- 
ical point of $ 4 theory. This rather bold conjecture led us 
to some surprising yet testable predictions for the scaling di- 
mensions of various operators in QED3. Future work, such as 
lattice simulations of QED3 14 and a higher-order renormaliza- 
tion group analysis assessing the stability of the QED3 fixed 
point, may provide further insight into the validity of this cor- 
respondence. 

To put our result in perspective, we have established du- 
ality between the two interacting conformal field theories by 
considering microscopic models, transforming to appropriate 
topological disorder variables, and studying detailed corre- 
spondence of the phases and transitions in the two systems. 
Our approach may also be useful in other contexts where such 
dualities have been conjectured42ii^ 

Another physically interesting setting in which this dual 
approach can be applied is the easy-plane spin- 1/2 triangu- 
lar antiferromagnet. In this case a direct Landau approach is 
not accessible, while the techniques introduced here permit an 
analysis in the dual language. Such an analysis can provide in- 
sight into critical spin-liquid phases, and will be carried out in 
a forthcoming paper. This may be useful for understanding re- 
cently observed spin liquids in the quasi-2D spin- 1/2 triangu- 
lar antiferromagnets 1 ' 3 CS2CUCI4 and k-(ET) 2 Cu2(CN)3. The 
easy-plane kagome antiferromagnet is another example of a 
frustrated spin system that might be fruitfully explored within 
a dual fermionic formulation. More generally, we expect that 
fermionization can provide an effective tool for describing 



APPENDIX A: APPLICATION OF FERMIONIZED 
VORTICES TO BOSONS AT v = 1 

In this Appendix, we provide some details of the fermion- 
ized vortex approach to bosons at filling factor v = 1 dis- 
cussed in the Introduction. This proceeds as for the easy-plane 
spins, which in fact we viewed as bosons in a magnetic field. 
To be closer to the continuum, we now imagine bosons say 
on a square lattice at low average density p <C 1 per site. 
The dual theory in the present v = 1 case has bosonic vor- 
tices at filling factor fduai = 1/V = L Vortices are at finite 
density because of the original magnetic field, while the dual 
magnetic field arises because the original bosons are at finite 
density. However, unlike the original bosons, vortices inter- 
act via a 2D electromagnetic interaction, and fermionization 
of vortices is more controlled. Upon Chern-Simons flux at- 
tachment and taking a flux-smeared mean-field state in which 
the statistical flux cancels the dual magnetic field, we obtain 
fermionic vortices in zero average field coupled to U(l) gauge 
field and Chern-Simons gauge field. This theory is very sim- 
ilar to Eq. Q, except that vortices are nonrelativistic at fi- 
nite density p, which is the same as the density of the original 
bosons. The fermionic vortex Lagrangian reads 



[V-i(q + A)] 7 

2777-vort 



i(a a + A )(^ip-p) 



—A • V x A A cxt ■ V x a 

Air 2tt 



(Al) 



Here m vor t is a vortex mass which we treat as a phenomeno- 
logical parameter. To allow computation of physical response 
properties, we have included an external probing field A ext 
minimally coupled to the original boson three-current Sj = 
Xr^. We also need to add the following bare action for the 
dual gauge field a: 



m 



Share M - 2(2 ^ )2 



J dr dr (Vao — d T a) 2 
^—tt [ dr dr dr' v rr > (V A a) r ( V A a) r 

>7T) Z J 



2(2tt) 

Here mo is roughly the original boson mass and v rT > encodes 
the boson density-density interaction. The limit ttiq — > cor- 
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responds to restricting the original bosons to the lowest Lan- 
dau level. With the above ingredients and using an RPA ap- 
proximation for integrating out the fermions, we can, e.g., ob- 
tain the physical response properties of the compressible state 
in agreement with Ref . [KJ 

From the above action, we see that the statistical gauge field 
A is 'screened' by the dual gauge field a. More formally, 
changing to a = a + A and integrating over the field A, we 
obtain 



2m v 



which is the equivalent of Eq. (|4) for nonrelativistic vortices. 
For simplicity, we do not write the quadratic terms for the 
gauge field a that arise from Share [a] ■ For finite mo and short- 
range v rr / these have the general structure shown in Eq. (0}. 

We now specialize to the lowest Landau level limit mo — ► 
0. The final action then describes fermions coupled to the 
gauge field a with no Chern-Simons term, which is similar to 
the result of Ref.lla 

Observe that in the present formulation the fermions are 
neutral by construction and do not couple directly to the prob- 
ing field A cxt . However, as the following calculation shows, 
fermions carry electrical dipole moments oriented perpendic- 
ular to their momentum k and of strength proportional to \k\. 
Indeed, from the action Eq. JA2t we obtain the following 
equations of motion, 



a tp'^jj 



P ■ 



(A3) 
(A4) 



(Incidentally, the first equation simply means that in the 
mo — > limit the field ao is completely soft and its fluctu- 
ations fix the vortex density; clearly, the vortex fermioniza- 
tion procedure is the simplest here.) Now, the original boson 
density is given by 

VAa VAA VAa . , + , _, , A „ 
<Ws = — = = -= + (i/>H - p) , (A5) 

Z7T Z7T Z7T 

and using the equations of motion we find 

i r d 2 k 



Spbos(q) 



•9/2 



(A6) 



Writing 5phos(q) = —iq • P(q)> we interpret this as a charge 
density created by a medium with bound charges and polar- 
ization 

As claimed, fermions carry dipole moment ~ Afe. Thus, the 
fermionized-vortex formulation reproduces the dipole picture 
of the v = 1 state of bosons. 



APPENDIX B: TRANSFORMATION PROPERTIES OF THE 
NEGATIVE-ENERGY DIRAC SEA 

In this Appendix we obtain the transformation properties 
given in Eq. J77I of the filled negative-energy Dirac sea in the 



presence of a charge q = ±1 monopole, \DS, q). We first note 
that the quoted action of 7f cml and C represents a convention 
for the phases of \DS, q). Time reversal sends the negative- 
energy Dirac sea of a monopole into the negative-energy Dirac 
see of an antimonopole, while particle-hole additionionally 
fills the two zero-energy modes of the antimonopole. One 
can thus use the specified transformations of \DS, +1) to 
define the phases of \DS, ±1), while the transformations of 



\DS,-1) follow from % 



ferm 



1 andC 2 = 1. 



We can also give a more explicit construction by consider- 
ing some general properties of the single-particle wavefunc- 
tions in the presence of a vector potential a xx , that gives rise 
to 2irq flux. The fermion hopping Hamiltonian is 



H t , q = -t v [e^x'e-Kx* +h.c. 

(XX') 



(Bl) 



This can be formally diagonalized by expanding d x in the en- 
ergy basis: 



E 



(B2) 



where a — 1,2 is the sublattice index, </>^(x) is a two- 
component wavefunction, and CE q annihilates the state with 
energy E in the presence of 2nq flux. The above sum is over 
all energy eigenstates. 

With our gauge choice, we have a^ x , = — a xx ,. Replacing 

dx,a —> ^aB^L a therefore sends H t . q — > H t .- q . This implies 
a relation between charge q and — q wavefunctions, 



<^(x) = ^*(x). (B3) 

Moreover, the transformation <i X Q , — > a^^d Xt /3 sends H t <q — > 
—Ht,q, which implies the following relation between wave- 
functions with energy E and — E, 



C E (x)=a*<(x). 



(B4) 



These properties allow us to determine how the operators 
CEq, and in turn \DS, q), transform under the internal sym- 
metries 7f crm and C. Consider fermionic time reversal first. 
According to Eq. ( I48> . 7f crm sends d x — > d x and q — > —q 
since the flux changes sign. We find using Eqs. (IB2> and (IB3> 
that 7f crm transforms cs.q — * CE- q - According to Table IIIII 
particle-hole sends d xa — * u^gd^g and also changes the sign 
of the flux. Using Eqs. JB2t through ( IB4> . we then find that 



ce,, 



cLe -q under C. 



Fermionic time reversal. Define \vac, q) to be the state with 
2iiq flux and all fermion modes unoccupied. Since the flux 
changes sign under time reversal, one can choose the phases 
of the fermion vacua such that T{ cim \vac, q) = \vac, —q). We 
write the filled negative-energy sea as 



\DS,q) 



£<0 



,t 

'Eq 



vac, q), 



(B5) 



where 7 6 K can be selected arbitrarily. Putting together the 
above results, we find that 



T lerm \DS,q) = \DS,-q). 



(B6) 
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Particle-hole. Under particle-hole, we expect physically 
that the vacuum with 2irq flux should transform to a state with 
— 2irq flux and all energy states filled. Using C 2 = 1, one can 
write without loss of generality 

C\vac, q) = e^ s J] J Et _ q \vac, -q) (B7) 

E 

by appropriately ordering the creation operators on the right 
side. Note that 6 is not arbitrary since the phases of \vac, q) 
were already fixed. We now choose the phase 7 so that the 
negative-energy sea transforms under particle-hole as follows, 

C\DS,q)=fl_ q fl_ q \DS,-q). (B8) 

We are unable to determine how the operators CEq trans- 
form under the lattice symmetries without knowing the ex- 
plicit forms of the wavefunctions. Nevertheless, under reason- 
able assumptions we can use various identities to determine 
how the negative-energy Dirac sea transforms under T$ T , and 
(up to overall signs) 1Z X and R„/ 3 . 

Translations. To deduce how the Dirac sea transforms un- 
der translations, we assume that \DS, q) is an eigenstate of the 
translation and rotation operators, 

T Sl \DS,q) = e^- 5r \DS,q) (B9) 
R«,z\DS,q) = e im «"/ 3 \DS,q). (BIO) 

Here, p g and m q are the total momentum and angular mo- 
mentum, respectively, of the Dirac sea. Define triangular lat- 
tice vectors Sri = x and Sr2 = — + ^y. The following 
operator relations are expected to hold, 

TsnRn/3 = ^7r/3^ 2 (Bll) 
Tsr 2 R n /3 = R^/^Ts ri Ts r2 (B12) 

since the left and right hand sides transform the lattice identi- 
cally. Requiring these relations to hold on the negative-energy 
Dirac sea, we find that \DS, q) must carry no overall momen- 
tum, i.e., 

T Sr \DS,q) = \DS,q). (B13) 



Modified reflections. The modified reflection 1Z X is a simple 
coordinate reflection composed with time reversal and does 
not change the direction of the flux. We expect the negative- 
energy Dirac sea to be an eigenstate of this operator. Further- 
more, commutation with particle-hole and time reversal can 
be used to show that the q = ±1 states have the same eigen- 
value ( x , which can either be +1 or —1: 

n x \DS,q) = C x \DS,q). (B14) 

To determine the sign £ x , one would need to know the wave- 
functions cj)q in Eq. dB2i to deduce how the operators CEq 

transform under 1Z X . 

Rotations. To obtain the angular momentum m q defined in 
Eq. ( IB 10t . we first note that commutation with time reversal 
immediately implies that m_ g = — m q . We restrict m q fur- 
ther using the observation that the composite operation 

O = R-} 3 CR* /3 C (B15) 

sends c? X Q — > ~d x a and leaves the flux unchanged. This 
can be readily verified from Table [H]] recalling that R n /3 
interchanges the two sublattices. Equation J B 2b then im- 
plies that O transforms CE q — > — Ce 9 . The filled Dirac sea 
is therefore an eigenvector of this operation with eigenvalue 

{—l) N f , where Njr is the total number of fermions resid- 
ing in the filled negative-energy sea. Assuming that gauge- 
invariant states have even total fermion number and noting 
that such states require one of the two zero-modes to be occu- 
pied, we conclude that in this case is odd. We thus have 
0\DS, q) =-\DS,q). This implies that 

R n/3 \DS,q) = ^/ze-^'^DS^q), (B16) 

where Q^/z — ±1 is an overall sign we are unable to determine 
without knowing the explicit wavefunctions <f>® . 
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